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I.  Introduction 


1.1  Background 

Over  the  years  the  DAPPA  research  program  to  develop  regional  discrimi¬ 
nation  procedures  has  produced  numerous  advances  in  our  understanding  of 
the  regions'  phase  signals  generated  by  seismic  sources.  However,  the  pro¬ 
gress  toward  determination  of  reliable  discriminants  has  been  slow  and  there 
have  been  no  true  breakthroughs  The  principal  hindrances  in  this  regard 
appear  to  be  lack  of  complete  theoretical  understanding  of  factors  affecting 
regional  phase  generation  and  propagation  and  inadequate  regional  phase 
data  from  Soviet  events  in  areas  of  interest  to  test  4he  portability  of  promising 
disenmination  procedures. 

Two  factors  suggest  that  it  may  now  be  possible  to  make  significant 
improvements  to  our  understanding  of  regional  discriminants  as  they  apply  to 
Soviet  events.  Principal  among  these  is  the  recent  availability  of  high-quality 
digital  seismic  data  from  stations  within  the  Soviet  Union  and  in  adjacent 
border  regions.  These  data  a^e  expected  to  be  particularly  important  for  moni- 
tonng  smaller  events  which  would  be  undetected  at  most  teleseismic  stations. 
The  second  .'actor  is  the  improvement  in  theoretical  studies.  Findings  from 
such  studies  are  expected  to  be  useful  in  identifying  new,  more-powerful 
regional  discrimination  measures  which  can  be  tested  on  relevant  event  data¬ 
bases  both  in  a  controlled  environment  like  the  NTS  region  and  in  the  region 
surrounding  the  Soviet  testing  areas  where  the  procedures  ultimately  need  to 
be  applied. 

1.2  Research  Summary 

Research  desebbed  in  this  report  summarizes  the  results  of  the  first  year’s 
effort  of  a  three  year  program.  This  program  is  designed  to  systematically 
evaluate  the  capability  of  regional  stations  to  d'seriminate  underground  nuclear 
explosions  from  earthquakes  and  non-nuclear  explosions  for  various  environ¬ 
ments  of  merest.  During  the  first  year  of  this  program,  efforts  have  focused  on 
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empirical  studies  utilizing  regional  data  from  the  Chinese  Digital  Seismic  Net¬ 
work  (CDSN)  in  the  People's  Republic  of  China  (PRC)  and  the  Incorporated 
Research  Institutions  for  Seismology  (IRIS)  stations  in  the  Soviet  Union. 
These  stations  provide  good,  strong  Pn,  Pg,  and  Lg  signals  from  underground 
nuclear  explosions  at  the  East  Kazakh  test  sites  and  earthquakes  in  the  region. 

In  previous  studies  we  had  found  that  triggered  seismic  stations  of  the 
CDSN  frequently  missed  the  Lg  signals  from  Soviet  underground  nuclear  tests. 
Only  the  nearest  station,  WMQ  at  Urumchi  (R  =  960  km),  appeared  to  have 
recorded  Lg  signals  consistently  from  the  Soviet  events.  However,  additional 
investigation  has  revealed  that  CDSN  station  HIA  (R  =  2910  km)  also  records 
Lg  signals  from  many  of  the  larger  East  Kazakh  explosions.  In  addition,  review 
of  the  data  recorded  at  the  Soviet  IRIS  stations  at  ranges  between  1380  km 
and  2880  km  reveals  that  they  also  record  strong  regional  signals,  in  some 
cases  to  magnitudes  as  low  as  4.0  mb  or  so.  We  have  found  that  the  regional 
signals  from  such  small  events  can  normally  be  enhanced  relative  to  noise  by 
band-pass  filtering  with  a  center  frequency  near  1  Hz. 

As  a  follow-up  to  preliminary  analyses  of  the  regional  signals  at  WMQ 
which  suggested  differences  in  the  relative  spectral  content  of  Lg  versus 
regional  P  for  explosions  and  earthquakes,  in  the  current  research  we  con¬ 
ducted  a  different  spectral  analysis  with  a  larger  sample  of  earthquake  and 
explosion  records.  The  analyzed  database  included  27  East  Kazakh  explo¬ 
sions  and  32  regional  earthquakes.  The  latter  were  selected  to  provide  a 
variety  of  propagation  paths  with  epicentral  distances  which  encompassed 
those  of  the  East  Kazakh  explosions.  These  studies  again  showed  that  Lg/P 
ratios  for  regional  earthquakes  were  relatively  enriched  at  high  frequencies 
compared  to  similar  explosions.  Comparison  of  events  with  similar  magnitudes 
indicated  that  regional  P  signals  appeared  to  have  similar  spectral  shapes 
while  the  Lg  signals  were  enriched  at  high  frequencies.  Lg/P  ratios  which  were 
intermingled  at  1  Hz  were  con  ;letely  separated  at  the  one  sigma  level  at  fre¬ 
quencies  of  2  Hz  and  greater  for  the  SR  explosions  and  regional  earthquakes. 
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Of  the  27  East  Kazakh  explosions  in  the  WMQ  database,  18  were  at 
Shagan  River  (SR)  and  9  were  at  Degelen  Mountain  (DM).  These  provided  an 
opportunity  to  study  differences  in  the  relative  signal  excitation  between  the  two 
source  areas.  We  found  from  our  analyses  that  the  regional  phase  signals 
appeared  to  show  different  effects  than  those  observed  by  Nuttli  (1987)  in  com¬ 
paring  Lg  amplitudes  with  teleseismic  P.  In  particular,  our  studies  suggest  that 
the  regional  P  signals  are  relatively  stronger  from  SR  explosions  than  from  DM 
explosions  with  similar  Lg  signals.  This  appears  to  be  the  opposite  of  Nuttli's 
finding  in  comparing  teleseismic  P  from  the  two  sites.  We  also  analyzed  the 
spectral  behavior  of  the  P  and  Lg  signals  recorded  at  WMQ  from  DM  versus 
SR  explosions.  These  results  revealed  a  tendency  for  the  DM  explosions  to 
have  Lg/P  spectral  ratios  which  were  slightly  more  earthquake-like  (i.e.  larger 
Lg/P  at  high  frequencies),  but  not  enough  for  the  source  site  differences  to 
jeopardize  the  discrimination  capability  of  the  measurement. 

An  issue  which  needs  further  study  is  the  degree  to  which  propagation 
path  differences  might  affect  the  Lg/P  spectral  ratios.  Since  co-located  earth¬ 
quake  and  explosion  sources  do  not  exist,  we  are  forced  to  compare  events 
which  have  somewhat  different  propagation  paths  to  the  available  regional 
monitoring  stations.  In  principle,  such  comparisons  would  be  most  valid  if  the 
signal  measurements  were  compensated  for  attenuation  differences.  To  obtain 
some  insight  into  these  effects,  we  have  made  some  preliminary  computations 
of  the  Q  values  for  Lg  signals  propagating  from  the  East  Kazakh  test  sites  to 
the  regional  CDSN  and  Soviet  IRIS  stations.  These  Q  values  were  derived 
from  an  analysis  of  the  Lg  magnitude  residuals  and  gave  values  ranging  from  a 
low  value  of  428  for  station  GAR  to  a  high  value  of  761  for  station  OBN. 

1.3  Report  Organization 

The  report  is  divided  into  six  sections  including  this  introduction.  Section  II 
summarizes  the  regional  network  used  in  these  studies  and  summarizes  the 
general  characteristics  of  the  regional  signals  recorded  at  the  CDSN  stations 
and  the  Soviet  IRIS  stations  from  East  Kazakh  explosions  and  regional  earth¬ 
quakes.  Section  HI  provides  information  on  the  spectral  analyses  performed  on 
the  WMQ  database  and  gives  the  results  of  the  Lg/P  spectral  ratio  measure¬ 
ments.  In  Section  IV  we  present  the  results  comparing  the  relative  Lg  and 
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regional  P  excitation  for  SR  and  DM  explosions.  Section  V  contains  the  prel¬ 
iminary  analysis  performed  on  the  Lg  magnitude  residuals  to  derive  effective  Q 
values  for  the  paths  to  the  regional  stations  from  the  East  Kazakh  test  sites. 
Finally,  Section  VI  provides  a  summary  and  identifies  future  research  plans. 
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II.  Database  for  CDSN  and  Soviet  IRIS  Stations 


2.1  The  Regional  Network 

As  noted  in  the  Introduction,  research  efforts  during  the  past  year  under 
this  contract  have  focused  on  analyses  of  the  regional  signals  recorded  at  sta¬ 
tions  in  the  CDSN  and  Soviet  IRIS  networks.  Figure  1  shows  the  locations  of 
these  stations  with  respect  to  the  principal  Soviet  underground  nuclear  test 
sites  at  Shagan  River  (SR)  and  Degelen  Mountain  (DM).  The  CDSN  stations 
range  in  distance  (cf.  Table  1)  from  about  960  km  for  station  WMQ  to  4380  km 
for  the  station  designated  QIO  on  the  map.  The  distance  to  the  test  sites  for 
the  Soviet  IRIS  stations  is  between  about  1380  km  for  station  GAR  and  2880 
km  for  station  OBN.  As  can  be  seen  from  the  map  in  Figure  1 ,  the  broadest 
gap  in  azimuthal  coverage  (=  125°)  is  to  the  north  of  the  test  site.  A  smaller 
gap  (=  80°)  occurs  to  the  south.  However,  it  should  be  noted  that  seismic  data 
are  not  currently  available  from  the  Soviet  IRIS  station  at  IRK  or  from  three  of 
the  CDSN  stations  designated  as  ENS,  QIO,  and  SSH  on  the  map;  so  actual 
regional  coverage  is  not  as  good  as  it  may  appear  from  the  station  map. 

The  stations  of  both  the  CDSN  and  Soviet  IRIS  networks  include  high- 
quality,  broad-band  three-component  systems  with  digital  recording  at  sampling 
rates  of  20  samples  per  second  and,  for  CDSN,  also  40  samples  per  second. 
The  stations  in  both  networks  also  include  low-frequency  elements  with  lower 
sampling  rates,  but  these  latter  have  not  been  analyzed  in  the  current  investi¬ 
gation.  An  important  difference  between  the  CDSN  and  Soviet  IRIS  recording 
systems  is  that  the  CDSN  is  a  triggered  system  while  the  IRIS  stations  record 
continuously.  As  a  result,  the  CDSN  stations  record  only  those  events  which 
produce  ground  motion  at  the  individual  stations  above  some  signal-to-noise 
threshold.  This  triggering  system  has  been  a  persistent  problem  for  recording 
longer-range  regional  events  since  this  network  first  became  operational.  In 
particular,  the  triggering  algorithm  is  not  optimized  for  detection  of  Lg  signals 
from  events  at  ranges  comparable  to  SR/DM.  Regional  signals  from  SR/DM 
explosions  below  about  4.8  mb  are  missed  at  even  the  nearest  CDSN  station 
(viz  WMQ);  and  at  the  distant  CDSN  stations  the  trigger  algorithm  usually 
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Figure  1.  Locations  of  CDSN  and  Soviet  IRIS  stations  with  re¬ 
spect  to  principal  East  Kazakh  test  sites  at  SR/DM. 
Solid  triangles  indicate  stations  from  which  useful 
Lg  signals  are  obtained. 
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Table  1.  Distance  From  Soviet  Test  Sites  For 
CDSN  and  Soviet  IRIS  Stations 


Station  Designation 

Distance  (km) 

Azimuth  (°) 

WMQ  (CDSN) 

960 

133 

GAR  (IRIS) 

1380 

212 

ARU  (IRIS) 

1530 

306 

IRK*(IRIS) 

1780 

72 

LZH  (CDSN) 

2530 

118 

KIV  (IRIS) 

2800 

270 

OBN  (IRIS) 

2880 

298 

HIA  (CDSN) 

2910 

76 

BJI  (CDSN) 

3100 

97 

ENS’(CDSN) 

3350 

120 

KMI  (CDSN) 

3440 

135 

MDJ  (CDSN) 

3800 

79 

SSH*(CDSN) 

4100 

105 

QIO*(CDSN) 

4380 

131 

No  useful  explosion  data  is  available  from  these  stations  at  this  time 


shuts-off,  for  even  the  largest  SR/DM  explosions,  prior  to  the  Lg  signal  and,  in 
most  cases,  doesn’t  reactivate.  Therefore,  the  analysis  of  the  CDSN  data 
presented  here  focuses  on  two  stations,  WMQ  and  HIA,  which  have  frequently 
recorded  Lg  signals  from  SR/DM  explosions.  The  other  CDSN  stations  nor¬ 
mally  record  only  the  regional  P  phases  from  SR/DM  explosions  and  occasion¬ 
ally  Lg  phases  from  earthquakes  nearer  individual  stations.  Analyses  of  the 
latter  data  will  be  included  in  subsequent  years  of  this  research  effort. 

With  regard  to  the  Soviet  IRIS  stations,  although  the  systems  are  designed 
to  record  continuously,  in  reality  this  goal  has  not  always  been  met  during  the 
initial  year  or  so  of  network  operation.  Station  down-time  has  been  a  recurring 
problem  particularly  at  some  of  the  stations.  As  a  result,  several  SR/DM  explo¬ 
sions  and  regional  earthquakes  have  been  missed  at  various  IRIS  stations. 
Furthermore,  during  this  start-up  phase  of  network  operation,  instrumentation 
and  gain  levels  have  occasionally  changed  and,  in  some  cases,  have  been 
incorrectly  reported.  This  has  provided  some  obstacles  to  the  determination  of 
absolute  ground-motion  levels  associated  with  the  regional  signals  recorded  at 
the  Soviet  IRIS  stations,  but  most  of  these  have  now  been  resolved. 

2.2  The  CDSN  Database 

The  high-quality  CDSN  network  has  been  in  routine,  nearly  continuous 
operation  since  late  1986.  Over  this  period  the  network  stations  have  recorded 
numerous  underground  nuclear  explosions  from  the  SR/DM  test  site.  For 
analyzing  the  characteristics  of  regional  signals  from  these  events,  we  have 
found  the  broadband  channels,  with  a  digitizing  rate  of  20  samples  per  second, 
to  be  most  useful.  This  band  appears  adequate  to  recover  the  spectral  charac¬ 
teristics  over  most  frequencies  for  which  the  signal-to-noise  level  is  above  1 .0 
for  the  SR/DM  explosions.  The  broadband  database  is  also  somewhat  more 
complete  than  the  short-period  band  (digitized  at  40  samples  per  second), 
which  experienced  more  of  the  triggering  problems  described  above. 

Figure  2  shows  the  broadband,  vertical-component  records  at  station 
WMQ  (R  =  960  km)  from  27  underground  nuclear  explosions.  Information  on 
the  locations  and  magnitudes  of  these  explosions  obtained  from  NEIS  are 
shown  in  Table  2.  Of  these  27  events  18  were  located  at  SR  and  9  at  DM. 
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Table  2.  Information  on  SR  &  DM  Explosions 
Recorded  at  CDSN  Stations 


Date 

Origin  Time 

Lat(°N) 

Lon(°E) 

Site 

m^j 

mb(Lg)  -  NOR 

WMQ 

HIA 

03/1 2/87 

01:57:17.2 

49.94 

78.82 

SR 

5.5 

5.22 

X 

_ 

04/03/87 

01:17:08.0 

49.93 

78.83 

SR 

6.2 

6.05 

X 

X 

05/06/87 

04:02:05.6 

49.83 

78.13 

DM 

5.6 

- 

X 

- 

06/06/87 

02:37:07.0 

49.86 

78.11 

DM 

5.3 

- 

X 

- 

06/20/87 

00:53:04.8 

49.91 

78.74 

SR 

6.1 

5.97 

X 

X 

07/17/87 

01:17:07.0 

49.80 

78.11 

DM 

5.8 

- 

X 

- 

08/02/87 

00:58:06.8 

49.88 

78.93 

SR 

5.9 

5.87 

X 

X 

11/15/87 

03:31:06.7 

49.87 

78.79 

SR 

6.0 

5.97 

X 

X 

12/13/87 

03:21:04.8 

49.99 

78.84 

SR 

6.1 

6.09 

X 

X 

12/20/87 

02:55:06.7 

49,83 

78.00 

DM 

4.8 

- 

X 

- 

12/27/87 

03:05:04.7 

49.86 

78.76 

SR 

6.1 

6.05 

X 

V 

02/06/88 

04:19:06.3 

49.80 

78.06 

DM 

5.0 

- 

X 

- 

02/13/88 

03:05:05.9 

49.95 

78.31 

SR 

6.1 

6.04 

X 

- 

04/03/88 

01 :33:05.8 

49.92 

78.95 

SR 

6.1 

6.06 

X 

X 

04/22/88 

09:30:06.9 

49.82 

78.12 

DM 

4.9 

- 

X 

- 

05/04/88 

00:57:06.8 

49.93 

78.77 

SR 

6.1 

6.04 

X 

- 

06/14/88 

02:27:06.4 

50.05 

79.01 

SR 

5.0 

- 

X 

- 

09/14/88 

04:00:00.0 

49.83 

78.81 

SR 

6.1 

5.97 

X 

X 

10/18/88 

03:40:06.6 

49.87 

78.08 

DM 

4.9 

- 

X 

- 

11/12/88 

03:30:03.7 

50.08 

78.99 

SR 

5.3 

5.26 

X 

- 

11/23/88 

03:57:06.7 

49.82 

78.07 

DM 

5.3 

- 

X 

- 

12/17/88 

04:18:06.9 

49.89 

78.93 

SR 

5.9 

5.80 

X 

X 

01/22/89 

03:57:06.6 

49.92 

78.83 

SR 

6.1 

5.96 

- 

X 

02/12/89 

04:15:06.8 

49.93 

78.74 

SR 

5.9 

5.79 

X 

X 

02/1 7/89 

04:01 :06.9 

49.87 

78.08 

DM 

5.0 

- 

X 

- 

07/08/89 

03:47:01.9 

50.66 

78.51 

SR 

5.6 

- 

X 

X 

09/02/89 

04:17:01.6 

50.85 

78.94 

SR 

5.4 

- 

X 

- 

10/19/89 

09:49:59.7 

50.47 

78.72 

SR 

6.0 

5.79 

X 

X 
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For  purposes  of  associating  these  events  with  a  particular  test  site,  we  have 
used  the  observation  that  the  event  locations  fall  into  two  clusters  with  the  DM 
cluster  located  about  50  km  west  of  the  SR  cluster.  The  events  cover  a  mag¬ 
nitude  range  from  4.8  to  6.2  mb  with  many  of  the  smaller  magnitude  explosions 
having  locations  at  DM.  East  Kazakh  explosions  with  magnitudes  below  about 
4.8  mb  apparently  do  not  trigger  the  detector  at  WMQ.  It  can  be  noted  in  Fig¬ 
ure  2  that  even  the  magnitude  4.8  mb  explosions  appear  to  have  large  signal- 
to-noise  levels;  and  it  should,  therefore,  be  possible  to  detect  regional  P  and  Lg 
signals  at  WMQ  from  much  smaller  East  Kazakh  explosions.  Although  some 
modifications  to  the  automatic  detection  algorithm  at  WMQ  have  been  made 
over  the  years  and  have  improved  acquisition  of  East  Kazakh  explosion  data,  it 
seems  clear  that  continuous  recording  or  changes  to  the  detector  could  enable 
recovery  of  much  smaller  events  at  WMQ.  One  or  two  larger  events  were  also 
not  recorded  at  WMQ  because  the  station  was  down  at  the  time  of  the  events. 

As  can  be  seen  in  Figure  2,  the  regional  P  and  Lg  signals  at  WMQ  from 
the  East  Kazakh  explosions  are  strong.  The  records  also  show  a  strong  Rg 
phase  whose  dispersion  characteristics  seem  to  be  remarkably  consistent  from 
event  to  event  but  whose  excitation  relative  to  the  other  regional  phases 
appears  to  be  somewhat  variable.  This  latter  observation  needs  further  investi¬ 
gation.  We  have  described  the  general  characteristics  of  the  regional  signals 
recorded  at  WMQ  in  a  previous  report  (cf.  Bennett  et  at.,  1989)  and  won't 
describe  them  further  here.  In  the  previous  report  we  also  described  spectral 
analyses  performed  on  the  regional  P  and  Lg  signals  at  WMQ  from  several 
East  Kazakh  explosions.  Spectral  analyses  of  the  Lg  signals  revealed  that  they 
were  peaked  in  a  relatively  narrow  frequency  band  between  0.3  and  1.C  Hz. 
The  regional  P  waves  showed  somewhat  broader  peaks  with  maxima  in  the 
band  from  about  0.5  to  3  Hz.  We  will  address  the  spectral  characteristics  of 
the  regional  signals  at  WMQ  more  fully  in  the  next  section  of  this  report. 

Our  earthquake  database  at  WMQ  currently  includes  32  events.  These 
are  listed  in  Table  3  and  their  locations  are  plotted  in  Figure  3.  The  earth¬ 
quakes  ranged  in  magnitude  from  4.3  to  5.9  mb.  Because  the  area  in  the 
immediate  vicinity  of  the  East  Kazakh  test  sites  is  generally  aseismic,  the  com¬ 
parative  earthquake  sample  was  selected  to  give  events  with  epicentral  dis¬ 
tances  which  bounded  the  explosion  distance  range.  In  addition,  some  events 
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Table  3.  Information  on  Earthquakes 
Recorded  at  WMQ 


Date 

Origin  Time 

Lat(°N) 

Lon(°E) 

mb 

R(km) 

11/18/86 

13:27:01.0 

40.06 

77.56 

4.7 

938 

12/14/86 

03:19:17.0 

47.31 

83.31 

5.0 

517 

12/20/86 

23:08:16.5 

36.75 

93.66 

5.3 

935 

01/05/87 

22:52:46.5 

41.96 

81.32 

5.9 

560 

01/24/87 

08:09:21.0 

41.53 

79.32 

5.9 

732 

01/24/87 

13:40:40.0 

41.44 

79.25 

5.2 

741 

03/05/87 

02:33:39.0 

35.41 

87.39 

4.5 

935 

04/09/87 

07:25:35.7 

35.50 

87.07 

4.8 

926 

04/09/87 

20:01:18.0 

35.59 

80.47 

4.9 

1103 

04/30/87 

05:17:37.0 

39.76 

74.57 

5.7 

1179 

06/08/87 

13:30:36.0 

39.79 

74.69 

5.1 

1169 

08/05/87 

10:24:21.0 

41.36 

82.11 

4.8 

534 

09/03/87 

09:08:12.0 

38.83 

75.32 

4.8 

1174 

09/16/87 

17:57:26.4 

52.09 

95.70 

4.8 

1096 

09/18/87 

21:58:41.0 

47.02 

89.66 

5.3 

387 

10/06/87 

13:06:20.3 

43.44 

88.55 

4.8 

81 

12/03/87 

23:51 :43.0 

39.54 

77.55 

4.7 

969 

12/06/87 

16:20:45.2 

37.44 

94.61 

4.7 

919 

12/17/87 

12:17:25.0 

41.94 

83.20 

5.1 

423 

12/22/87 

00:16:39.0 

41.36 

89.64 

5.9 

317 

01/02/88 

22:02:36.0 

40.06 

77.34 

4.9 

954 

01/09/88 

03:55:05.3 

39.09 

71.50 

5.4 

1449 

02/08/88 

17:49:19.8 

43.73 

83.76 

4.3 

317 

03/25/88 

02:07:56.0 

44.71 

79.60 

4.5 

654 

04/01/88 

01:27:16.2 

47.53 

89.65 

4.6 

439 

05/02/88 

02:13:26.0 

40.26 

82.20 

4.9 

603 

05/25/88 

00:05:23.0 

40.57 

77.62 

4.9 

907 

05/25/88 

18:21:58.0 

42.01 

85.69 

5.2 

259 

06/17/88 

13:30:45.0 

42.94 

77.50 

5.3 

832 

06/30/88 

15:25:15.5 

50.23 

91.14 

5.0 

759 

07/23/88 

07:38:09.0 

48.72 

90.51 

5.5 

586 

11/15/88 

16:56:45.9 

42.08 

89.40 

5.0 

238 
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Fiqure  3.  Locations  o  f  M  earthquakes  recorded  at  WM<.)  with  respect  to  So  vie 
nuclear  explosion  tests  at  SR/PM. 


at  smaller  and  larger  ranges  from  WMQ  were  selected  to  enable  a  more  com¬ 
plete  analysis  of  attenuation  and  other  path-dependent  effects  for  the  region. 
As  a  result  the  earthquakes  cover  an  epicentrai  distance  range  from  about  80 
km  to  1450  km.  The  earthquakes  recorded  at  WMQ  typically  show  strong  Lg 
signals  which,  in  most  cases,  are  larger  than  the  frequently  emergent  P 
phases.  Some  of  these  characteristics,  based  on  a  smaller  data  sample,  were 
described  in  our  previous  report  (cf.  Bennett  et  al.,  1989)  and  examples  of  the 
observed  signals  were  presented.  In  Section  III  of  this  report,  we  describe  the 
results  of  spectral  analysis  performed  on  the  regional  earthquake  signals 
recorded  at  WMQ  and  identify  diagnostic  differences  in  the  signal  behavior 
between  the  earthquake  and  explosion  samples. 

Figure  4  shows  the  broadband,  vertical-component  records  at  station  HIA 
(R  =  2910  km)  from  12  SR  underground  nuclear  explosions.  Information  on 
the  locations  and  magnitudes  of  these  explosions  is  presented  in  Table  2 
above.  The  events  recorded  at  HIA  cover  a  very  limited  range  of  magnitudes 
between  5.6  and  6.2  mb.  Smaller  magnitude  explosions  at  the  East  Kazakh 
test  sites  apparently  do  not  trigger  the  detector.  From  the  broadband  'ecord- 
ings  in  Figure  4,  it  appears  that  the  magnitude  5.6  event  has  a  signal-to-noise 
level  of  about  3:1.  Some  enhancement  of  the  recording  capability  at  HIA  for 
East  Kazakh  explosions  again  might  be  achieved  by  continuous  recording  or 
modifying  the  triggering  algorithm,  but  the  signal-to-noise  level  suggests  that  it 
may  not  be  possible  to  push  the  threshold  at  HIA  below  about  magnitude  5.0 
mb.  It  should  also  be  noted  in  Figure  4  that  the  records  of  earlier  explosions 
obtained  from  HIA  consist  of  interrupted  P  and  Lg  segments.  The  intervening 
trace  dead  time  corresponds  to  a  time  interval  during  which  the  automatic 
detector  shuts  off  before  retriggering  on  the  Lg  arrival.  This  problem  appears 
to  have  been  corrected  for  the  more  recent  explosions  producing  the  continu¬ 
ous  records  shown  for  the  bottom  five  events. 

The  SR  explosion  records  obtained  at  HIA  in  Figure  4  show  relatively 
strong  Lg  signals  following  a  relatively  weak  P  phase.  The  relative  phase 
characteristics  are  observed  to  be  very  consistent  from  event  to  event.  The  Lg 
window  includes  a  higher-frequency  segment  followed  by  a  dispersed  packet  of 
low-frequency  energy  apparent  on  most  records.  Typical  P  and  Lg  spectra  for 
SR  explosions  at  HIA  are  shown  in  Figure  5.  As  with  the  WMQ  signals,  the  Lg 


14 


4199  *T  WIH 

•  IZ2»T  41/C  /» 


5 


Fiqure  4.  Vertical-component  records  at  HTA  from  12  SR  explosions 


spectra  at  HIA  are  sharply  peaked  at  a  frequency  less  than  1  Hz  (viz  between 
about  0.3  and  0.8  Hz).  There  is  some  evidence  of  a  secondary  peak  in  the  Lg 
spectra  at  a  somewhat  higher  frequency  (around  4  Hz);  but  the  signal-to-noise 
level  in  this  band  is  low  and  the  peak  could  be  related  to  noise.  The  Lg  spec¬ 
tra  are  generally  at  or  near  the  noise  level  for  frequencies  of  2  Hz  and  higher. 
The  regional  P-wave  spectra  from  the  East  Kazakh  explosions  recorded  at  HIA 
are  peaked  at  a  higher  frequency  (around  1  Hz)  than  the  Lg.  Furthermore,  the 
P-wave  spectra  are  somewhat  broader  than  the  Lg  with  relatively  more  energy 
at  frequencies  above  1  Hz.  One  peculiarity  of  the  Lg  signal  spectrum  for  the 
September  14,  1988  explosion  at  SR  is  its  failure  to  return  to  the  noise  level  at 
high  frequencies.  We  are  attempting  to  ascertain  whether  this  is  somehow 
related  to  the  signal  processing,  but  most  obvious  causes  have  been  elim¬ 
inated. 

2.3  The  Soviet  IRIS  Database 

The  network  of  Soviet  IRIS  stations  began  operation  in  September,  1988. 
As  noted  above,  the  operation  of  individual  stations  in  the  network  has  been 
somewhat  sporadic  during  the  initial  start-up  phase.  This  has  resulted  in  a 
failure  to  record  several  explosions  at  the  East  Kazakh  test  sites  because  sta¬ 
tions  were  not  operating  at  these  times.  Nevertheless,  excellent  regional 
seismic  records  have  been  obtained  from  several  Soviet  IRIS  stations  starting 
with  the  Soviet  JVE  shot  on  September  14,  1988.  Figure  6  shows  the  avail¬ 
able  vertical-component  records  obtained  from  nine  explosions  through  July, 
1989.  Two  additional  explosions  since  that  time  are  in  the  process  of  being 
installed  into  the  CSS  database.  Table  4  shows  the  location  and  magnitude 
information  on  the  nine  explosions.  Seven  of  the  events  were  located  at  SR 
and  two  at  DM.  The  magnitudes  range  from  3.8  to  6.1  mb.  The  smallest  event 
magnitude  (viz  3.8  mb)  was  estimated  by  Hansen  et  at.  (1990)  based  on  NOR- 
SAR  mb  measurements  with  a  bias  correction  applied  to  adjust  to  a  large, 
global-network  mb.  All  nine  of  the  explosions  were  recorded  at  ARU,  and  only 
one  other  Soviet  IRIS  station  (viz  GAR)  recorded  more  than  half  of  the  events. 
This  is  not  a  detection  problem;  the  missing  stations  were  simply  not  operating. 

As  can  be  seen  in  Figure  6,  the  Soviet  IRIS  stations  generally  record 
strong,  rather  complex  regional  P  phases  and  equally  strong  Lg  on  these 
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Table  4.  Information  on  SR  &  DM  Explosions 
Recorded  at  Soviet  IRIS  Stations 


Date 

Origin  Time 

Lat(°N) 

Lon(°E) 

Site 

m*) 

mB(lg)  -  NOR 

GAR 

ARU 

KIV 

OBN 

09/14/88 

04:00:00.0 

49.83 

78.81 

SR 

6.1 

5.97 

X 

X 

X 

x  ! 

09/26/88* 

07:45:00.0 

49.90 

78.00 

SR 

3.8 

- 

- 

X 

- 

10/18/88 

03:40:06.6 

49.87 

78.08 

DM 

4.9 

- 

X 

X 

- 

- 

11/12/88 

03:30:03.7 

50.08 

78.99 

SR 

5.3 

5.26 

X 

X 

- 

- 

12/17/88 

04:18:06.9 

49.89 

78.93 

SR 

5.9 

5.80 

- 

X 

- 

X 

01/22/89 

03:57:06.6 

49.92 

78.83 

SR 

6.1 

5.96 

- 

X 

- 

- 

02/12/89 

04:15.06.8 

49.93 

78.74 

SR 

59 

5.79 

- 

X 

- 

- 

02/17/89 

04.01:06.9 

49.87 

78.08 

DM 

5.0 

- 

X 

X 

- 

07/08/89 

03:47:01.9 

50.66 

78.51 

SR 

5.6 

- 

X 

X 

X 

X 

‘Source  intormation  on  (his  event  is  from  Hansen  et  al.  (1990);  the  location  and  origin  time  are  approxi¬ 
mate. 


19 


broadband  systems.  The  only  exception  with  regard  to  Lg  is  station  KIV  for 
which  the  propagation  path  from  the  East  Kazakh  test  sites  crosses  the  north¬ 
ern  Caspian  Sea,  which  has  previously  been  recognized  as  blocking  Lg  propa¬ 
gation  ^cf.  Piwinskn  and  Spnnger,  1978,  Bennett  ot  si.,  1981).  The  Lg  signal  at 
KIV  in  the  records  of  Figure  6  appears  to  be  at  the  P-coda  level  and  indistinct. 
In  contrast,  Lg  signal  levels  at  the  two  nearest  IRIS  stations  (ARU  and  GAR) 
from  East  Kazakh  explosions  are  particularly  strong  for  even  small  magnitude 
events.  Although  this  observation  may  not  be  immediately  obvious  from  the 
broadband  records  in  Figure  6,  we  will  show  later  that  even  the  five  records 
from  the  smaller  magnitude  events  at  the  bottom  of  the  figure  have  signal-to- 
noise  levels  well  above  1 .0  over  some  frequency  bands. 

Spectral  analyses  were  performed  on  waveform  segments  corresponding 
to  P,  Lg,  and  pre-P  noise  for  each  of  the  records  in  Figure  6.  Examples  of  the 
spectra  for  two  explosions  at  each  station  are  shown  in  Figures  7  and  8.  The 
analyses  revealed  that  the  energy  in  the  Lg  spectra  is  concentrated  in  a  narrow 
frequency  band  which  peaks  at  a  frequency  less  than  1.0  Hz.  The  P-wave 
spectra  at  the  various  Soviet  IRIS  stations  are  generally  much  broader  and 
remain  above  the  noise  level  to  higher  frequencies.  The  noise  spectra  them¬ 
selves  show  a  null  at  a  frequency  near  1  Hz.  This  seems  to  be  fairly  con¬ 
sistent  between  all  stations  and  events.  The  noise  spectra  reach  a  peak  at  a 
somewhat  higher  frequency  before  declining  again  above  5  Hz,  as  the 
antialiasing  filter  in  the  IRIS  recording  system  takes  effect.  The  spectral 
behavior  of  the  P  and  Lg  signals  at  the  Soviet  IRIS  stations  is  in  general  con¬ 
sistent  with  what  would  be  expected  from  an  attenuation  model;  spectral  decay 
with  frequency  above  the  peak  values  is  more  rapid  at  the  more  distant  sta¬ 
tions  because  of  higher  attenuation.  As  a  result  of  this  attenuation,  the  Lg  sig¬ 
nals  in  particular  are  dominated  by  low  frequencies  at  the  more-distant  sta¬ 
tions.  It  can  be  seen  from  the  spectral  plots  in  Figure  8  that  the  Lg  signals  at 
OBN  and  KIV  have  low  signal-to-noise  levels  above  about  2  Hz  for  even  large 
explosions  at  the  East  Kazakh  test  sites.  It  is  apparent  from  this  observation 
that  to  resolve  high-frequency  characteristics  of  regional  Lg  signals  from  most 
Soviet  underground  nuclear  explosions  will  require  regional  stations  at  ranges 
much  closer  than  those  of  OBN  and  KIV  (i.e.  R  «  2900  km).  The  regional  P- 
wave  spectra  in  Figure  8  indicate  a  somewhat  broader  band  of  useful  frequen¬ 
cies.  For  the  P-wave  windows  signal-to-noise  levels  are  greater  than  1.0  to 
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Figure  7.  Regional  P,  Lg  and  noise  spectra  measured  at  Soviet 
IRIS  stations  from  the  SR  explosion  of  09/14/88 
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frequencies  of  5  Hz  or  more  at  the  more  distant  stations  and  to  1 0  Hz  at  the 
two  nearer  stations. 

One  clear  result  of  the  spectral  analyses  of  the  regional  signals  from  East 
Kazakh  explosions  measured  at  Soviet  IRIS  stations  is  that  the  signal-to-noise 
level  is  maximum  at  a  frequency  of  about  1  Hz  for  both  Lg  and  regional  P. 
Although  this  frequency  range  may  vary  slightly  from  station  to  station  and  may 
be  somewhat  less  for  Lg  (cf.  Figure  8),  the  observation  suggests  that  useful 
regional  signals  can  be  extracted  from  the  background  noise  by  bandpass 
filtering  at  a  frequency  of  about  1  Hz.  This  is  illustrated  in  Figure  9  where  we 
have  applied  a  narrow  band-pass  filter  with  corners  at  0.8  and  1 .6  Hz  to  the 
ground  motion  time  histories  recorded  at  station  ARU  for  several  underground 
nuclear  tests  at  the  East  Kazakh  test  sites.  The  events  shown  ranged  in  mag¬ 
nitude  from  6.1  mb  to  3.8  mb.  The  raw,  unfiltered  traces  are  shown  on  the  left 
and  the  filtered  traces  are  on  the  right.  It  can  be  seen  that,  for  the  magnitude 
6.1  and  5.9  explosions,  the  signal-to-noise  level  is  good  prior  to  application  of 
the  filter.  For  these  events  the  filter  reduces  some  of  the  longer-  period  energy 
arriving  late  in  the  Lg  or  at  the  beginning  of  the  Rg  windows.  For  the  magni¬ 
tude  5.3  mb  explosion,  the  signal  is  visible  in  the  raw,  broad-band  record  in 
spite  of  rather  strong  lower-frequency  noise.  The  signal  appearance  is  greatly 
enhanced  by  the  filtering  operation.  The  signal-to-noise  improvement  is  even 
more  dramatic  for  the  three  smaller-magnitude  events.  The  raw  traces  for  the 
magnitude  5.0  and  4.9  mb  explosions  show  little  evidence  of  signal  except  for  a 
few  higher-frequency  segments.  When  filtered  these  high-frequency  segments 
are  greatly  enhanced  and  revealed  to  be  dear  regional  P,  S,  and  Lg  signals. 
The  smallest  magnitude  explosion  for  which  Soviet  IRIS  data  are  available  was 
estimated  to  have  a  magnitude  of  3.8  mb  by  Hansen  et  al.  (1990).  Though  the 
magnitude  of  this  event  is  somewhat  questionable,  since  it  was  based  on  a  sin¬ 
gle  amplitude  measurement,  it  seems  unlikely  from  the  relative  amplitude  of  the 
Lg  that  this  explosion  had  a  magnitude  much  greater  than  4  mb.  For  this  mag¬ 
nitude  3.8  event,  the  raw  trace,  on  the  left  in  Figure  9,  shows  very  little  evi¬ 
dence  of  signal.  The  filtered  trace  on  the  right,  however,  shows  Lg  signals  with 
amplitudes  more  than  three  times  the  noise  level.  The  regional  P  signal  for 
this  explosion,  however,  is  not  much  bigger  than  the  ambient  noise  even  after 
filtering.  We  performed  similar  analyses  on  the  records  at  GAR  for  the  lower- 
magnitude  explosions  in  Figure  7  and  found  that  they  also  showed  clear 
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regional  signals  when  bandpass  filtered  at  about  1  Hz. 

During  the  first  year  of  this  contract,  we  have  recovered  only  a  small 
amount  of  earthquake  data  recorded  at  the  Soviet  IRIS  stations.  Soviet  IRIS 
stations  GAR  and  KIV  are  located  in  active  seismic  zones  along  the  southern 
Soviet  border  region  and  have  recorded  numerous  local  events.  However,  in 
our  discrimination  efforts  we  have  attempted  to  focus  on  earthquakes  which 
are  similar  in  magnitude  or  have  comparable  paths  to  the  available  explosions 
in  the  database.  These  restrictions  severely  limit  the  available  data  and  may 
need  to  be  relaxed  to  increase  the  database  of  regional  signals  from  natural 
events  recorded  at  the  Soviet  IRIS  stations. 

Two  prominent  sources  of  regional  signals  from  larger  earthquakes,  with 
magnitudes  comparable  to  Soviet  underground  nuclear  explosions,  during  the 
relatively  recent  operating  history  of  the  IRIS  stations  have  been  the  12/07/88 
earthquake  in  Armenia  and  its  aftershocks  and  also  several  larger  earthquakes 
in  the  Pamir-Hindu  Kush  region.  We  have  recovered  Soviet  IRIS  station  data 
for  some  of  these  events.  Initial  indications  are  that  the  earthquakes  produce 
strong  Lg  signals  which  frequently  have  larger  amplitudes  than  the  regional  P 
phases  on  the  20  samples  per  second  channels  at  the  IRIS  stations.  In 
several  cases  the  dominant  frequencies  of  the  Lg  signals  from  these  earth¬ 
quakes  appear  to  occur  at  lower  frequencies  than  comparable  explosions. 
This  is  particularly  true  for  some  of  the  larger-magnitude  earthquakes.  It 
should  be  noted,  however,  that  the  propagation  paths  for  these  events  are  also 
different  than  for  the  explosions.  In  future  analyses  we  plan  to  take  into 
account  such  propagation  differences  and  make  appropriate  adjustments  for 
attenuation  effects  which  could  alter  the  spectral  content  of  the  regional  P  and 
Lg  signals.  In  these  analyses  we  hope  to  incorporate  a  much  larger  sample  of 
regional  earthquakes  recorded  at  all  the  available  Soviet  IRIS  stations  which 
will  permit  more  complete  investigation  of  path  effects. 
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III.  Analyses  of  Spectral  Differences  Between 
Explosions  and  Earthquakes  Recorded  at  WMQ 


3.1  Objectives  of  Discrimination  Analyses 

The  primary  goal  of  this  research  program  is  to  develop  reliable  regional 
discrimination  procedures  which  can  be  applied  to  identify  seismic  events  in 
the  Soviet  Union.  The  focus  of  our  studies  to  date  under  this  contract  has 
been  empirical  observations  and  analyses  of  the  characteristics  of  regional 
phase  signals  from  Soviet  underground  nuclear  explosions  and  comparable 
earthquakes  recorded  at  high-quality  digital  seismic  stations  in  the  regional  dis¬ 
tance  range.  Such  data  are  likely  to  be  extremely  important  to  the 
identification  of  small  events  under  a  low-threshold  testing  limit.  Over  the 
years  several  time-domain  and  spectral  measures  of  seismic  signals  have 
been  proposed  as  potential  regional  discriminants  (cf.  Blandford,  1981; 
Pomeroy  et  al.,  1982).  However,  these  regional  discrimination  techniques  have 
not  been  completely  reliable;  and,  furthermore,  capability  to  extend  the  tech¬ 
niques  to  aseismic  regions  of  the  Soviet  Union  (including  the  principal  Soviet 
test  sites  in  East  Kazakh)  has  not  been  demonstrated.  The  latter  situation  is 
partly  due  to  limitations  on  our  theoretical  understanding  of  regional  phase  gen¬ 
eration  (particularly  for  explosion  sources)  and  a  lack  of  high-quality  seismic 
data  at  regional  distances  from  events  in  the  Soviet  Union.  It  is  only  recently 
that  such  data  have  become  available  with  the  advent  of  the  CDSN  stations  in 
China  and  the  even  more  recent  installation  of  the  Soviet  IRIS  stations.  In  a 
previous  report  (cf.  Bennett  et  al.,  1989)  we  presented  preliminary  results 
describing  general  characteristics  of  the  regional  signals  measured  at  a 
selected  CDSN  station,  WMQ,  from  explosions  at  the  East  Kazakh  test  sites 
and  nearby  earthquakes.  We  also  presented  a  band-pass  filter  analysis  of 
those  data  which  indicated  diagnostic  differences  in  the  relative  amplitudes  of 
the  Lg  and  regional  P  phases  at  certain  frequencies  from  explosions  and  earth¬ 
quakes  recorded  at  WMQ.  In  the  current  study  we  have  attempted  to  follow-up 
these  observations  with  more  thorough  spectral  analyses  which  were  per¬ 
formed  on  a  larger  data  sample. 
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3.2  Preliminary  Observations  of  Spectral  Behavior 


CDSN  station  WMQ  is  the  nearest  station  to  the  principal  Soviet  test  sites 
at  East  Kazakh  for  which  high-quality  digital  data  are  routinely  available.  As 
pointed  out  in  the  preceding  section  of  this  report,  WMQ  records  strong 
regional  P  and  Lg  signals  from  Soviet  underground  nuclear  tests  at  East 
Kazakh  (R  =  960  km)  and  also  from  earthquakes  at  epicentral  distances  which 
encompass  the  range  to  the  East  Kazakh  test  sites.  The  instrument  response 
for  the  WMQ  records  used  in  our  analyses  were  broadband,  and  we  applied 
different  narrow  band-pass  filters  to  derive  Smax/Pmax  and  Lg/Pmax  ratios  as  a 
function  of  frequency.  Figure  10,  taken  from  our  previous  report,  shows 
Sm^^max  ratios  measured  from  the  raw  broadband  vertical-component  records 
for  a  sample  consisting  of  12  East  Kazakh  explosions  and  20  nearby  earth¬ 
quakes.  The  earthquakes  were  basically  a  subset  of  those  shown  above  in 
Figure  3  and  Table  3  with  magnitudes  between  4.5  and  5.9  mb  and  ranges 
between  423  and  1176  km.  The  earthquakes  were  located  mainly  along  the 
southern  Soviet  border,  so  that  their  paths  to  WMQ  were  somewhat  different 
from  those  of  the  explosions.  This  raises  some  concern  about  the  effects  of 
dissimilar  paths  on  the  signal  measurements.  Nevertheless,  the  results  indi¬ 
cated  a  fairly  large  separation  by  about  a  factor  of  four,  on  average,  between 
explosions  and  earthquakes  with  the  explosion  ratio  being  much  smaller.  A 
few  anomalous  events  at  lower  magnitudes  had  Smax/Pmax  ratios  which  were 
intermingled  between  earthquakes  and  explosions. 

The  bandpass  filter  analyses  were  performed  to  explore  this  observation 
more  fully  and  possibly  to  increase  the  separation  between  the  two  source 
types.  Figure  11  shows  the  results  of  the  preliminary  analyses  which  were 
applied  to  10  East  Kazakh  explosions  and  13  earthquakes.  The  Lg/Pmax  ratios 
determined  from  the  filtered  data  for  explosions  and  earthquakes  are  observed 
to  be  intermingled  at  filter  center  frequencies  of  1  Hz  and  lower  but  separated 
in  all  cases  above  3  Hz.  The  average  Lg/Pmax  ratio  was  found  to  be  only  about 
30  percent  larger  for  earthquakes  than  for  explosions  in  the  low-frequency 
passband.  However,  the  explosion  ratio  fell  off  much  faster  toward  higher  fre¬ 
quencies  and,  as  a  result,  was  about  a  factor  of  ten  smaller  than  the  average 
earthquake  ratio  at  high  frequencies.  We  concluded  that  care  needed  to  be 
exercised  in  dealing  with  regional  discriminants  which  rely  on  signal  differences 
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in  limited  frequency  bands.  Some  recording  systems  may  have  responses 
which  are  too  narrow  to  resolve  the  differences.  For  example,  a  system  with  a 
relatively  narrow  response  peaked  near  1  Hz,  such  as  WWSSN,  might  not  see 
the  differences  in  the  Lg/Pmax  ratios  which  were  found.  This  might  be  an  expla¬ 
nation  why  Lg/Pmax  ratio  observations  have  not  always  provided  reliable 
discrimination  of  explosions  and  earthquakes. 

3.3  Further  Spectral  Analysis  of  WMQ  Data 

In  the  current  study  we  performed  an  alternate  spectral  analysis  on  the 
regional  P  and  Lg  signals  from  a  larger  sample  of  East  Kazakh  explosions  and 
earthquakes  recorded  at  WMQ.  Instead  of  narrow  band-pass  filtering,  discrete 
Fourier  transforms  were  computed  for  the  regional  P  and  Lg  signals;  and  spec¬ 
tral  ratios  were  derived  from  those.  For  these  analyses  we  used  the  vertical- 
component  signals  which  were  described  in  Section  II.  As  can  be  seen  in  Fig¬ 
ure  2  above,  the  regional  P  signals  from  East  Kazakh  explosions  recorded  at 
WMQ  are  quite  complex  with  apparent  multiple  phases  of  approximately 
equivalent  amplitudes  on  the  broadband  recordings.  The  earthquake  P 
phases,  although  normally  smaller  and  more  variable,  were  also  seen  to  be 
complex  in  most  distance  intervals  apparently  indicating  multiple  P  paths  in  the 
crust  and  upper  mantle. 

From  the  WMQ  data  identified  in  Tables  2  and  3  above  a  sample  of  19 
East  Kazakh  explosions  and  27  regional  earthquakes  were  analyzed.  The 
earthquakes  were  chosen  to  provide  a  distance  range  which  encompassed  the 
distances  from  WMQ  to  the  East  Kazakh  test  sites,  but  several  different 
azimuths  ere  represented  in  the  data  sample.  It  was  anticipated  that,  by  allow¬ 
ing  these  kinds  of  variations  in  propagation  path,  bias  associated  with  particular 
paths  could  be  avoided  or  at  least  diminished.  We  also  performed  some  prel¬ 
iminary  tests  on  the  P-wave  window  length  in  an  attempt  to  establish  a  stable 
measurement  scheme.  Pre-P  noise  spectra  were  also  calculated  for  each 
record  to  verify  the  frequency  range  over  which  the  individual  phase  spectra 
and  the  corresponding  ratios  could  be  considered  valid. 

The  initial  P-wave  arrivals  were  timed  using  a  phase  picker  based  on  the 
Stewart-Alien  algorithm  (cf.  Allen,  1982).  The  automatic  picks  were  then 
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verified  by  the  analyst.  The  data  were  of  sufficiently  high  quality  that  the 
automatic  picker  failed  to  find  the  correct  P-wave  arrival  time  for  only  one 
earthquake.  Time  marks  were  placed  on  the  records  corresponding  to  the  P- 
arrival  times  and  to  the  group  velocity  interval  between  3.5  and  2.9  km/sec, 
which  includes  most  of  the  Lg  energy  for  these  events.  The  pre-P  noise  time 
segment  was  defined  as  being  from  1  second  before  the  initial  P  arrival  to  30 
seconds  before  that  time.  In  cases  where  the  P  arrival  occurred  less  than  30 
seconds  after  the  record  start,  the  noise  window  was  taken  from  the  beginning 
of  the  record. 

In  order  to  avoid  spurious  high-frequencies,  each  window  was  tapered 
using  a  Hanning  window  applied  over  5  percent  of  each  end  of  the  time  seg¬ 
ment.  The  discrete  Fourier  transforms  were  then  calculated  for  each  window 
producing  amplitude  and  phase  spectra.  The  spectra  were  not  always  con¬ 
sistent  in  length  or  frequency  spacing  because  of  the  varying  sizes  of  the  time 
windows,  so  they  were  all  resampled  to  120  evenly  spaced  frequencies  using 
the  interpolation  method  of  Wiggins  (1976).  The  spectra  were  then  convolved 
with  a  boxcar  function  10  points  wide  to  provide  some  smoothing.  All  spectral 
ratios  were  calculated  using  these  resampled,  smoothed  spectra. 

3.4  Comparison  of  P-Wave  Windows 

As  pointed  out  above,  the  regional  P-wave  signals  recorded  on  the  broad¬ 
band  records  at  WMQ  from  East  Kazakh  explosions  and  regional  earthquakes 
are  frequently  complex.  The  vertical-component  records  in  Figure  2  above  for 
East  Kazakh  explosions  show  two  or  three  distinct  regional  P  phases  rising 
above  the  general  P  coda  level.  To  test  the  sensitivity  of  the  spectral  meas¬ 
urements  to  the  P-wave  window,  we  analyzed  the  P-wave  signals  for  two 
different  windows,  both  starting  1  second  prior  to  the  detected  P  arrival.  The 
short  P-wave  window  extended  to  10  seconds  after  the  P  arrival  and  the  long 
window  to  25  seconds  after  the  P  arrival.  The  short  window  included  mainly 
the  energy  in  the  first  distinct  P  phase  seen  in  the  explosion  records  of  Figure 
2  above.  The  long  P  window  appears  to  include  most  of  the  P-wave  energy 
seen  on  these  records. 
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Figures  12  and  13  show  the  analyses  performed  on  the  P-wave  signals  for 
two  East  Kazakh  explosions:  one  at  SR  (viz  03/i?/87)  and  one  at  DM  (viz 
05/06/87).  The  expanded  playouts  of  the  P  phases  at  the  tops  of  each  figure 
illustrate  the  complexity  of  the  regional  signals.  Several  arrivals  are  apparent 
in  the  first  25  seconds  of  the  signals.  One  interesting  feature  in  comparing  Fig¬ 
ures  12  and  13  is  that  the  initial  P  for  the  DM  event  appears  to  have  a  some¬ 
what  lower  predominant  frequency,  even  though  the  events  have  similar  body- 
wave  magnitudes  (viz  5.5  for  the  03/12/87  explosion  and  5.6  for  the  05/06/87 
explosion).  We  will  address  differences  in  the  regional  signals  between  SR 
and  DM  events  more  completely  in  a  subsequent  section  of  this  report. 

The  plots  in  the  lower  left  corners  of  Figures  12  and  13  show  the  spectra 
computed  for  these  events.  Spectra  are  shown  for  both  the  short  and  long 
regional  P  windows,  the  Lg  windows  and  the  noise.  In  addition,  the  lower  solid 
line  shows  the  ratio  of  the  spectra  for  the  short  P-wave  window  to  the  long  win¬ 
dow.  First,  it  can  be  noted  that  the  signal  spectra  generally  lie  well  above 
noise  over  the  entire  band  of  frequencies  plotted  here.  It  should  be  pointed  out 
that  the  instrument  response  has  not  been  removed  from  the  spectra 
presented  in  this  section.  Focusing  on  the  P-wave  spectra,  the  short-  and 
long-window  spectra  are  seen  to  match  quite  closely  at  high  frequencies;  but 
the  short-window  spectra  fall  below  the  long-window  spectra  at  frequencies 
less  than  about  3  Hz.  Thus,  the  general  tendency  is  for  the  long-window  spec¬ 
tra  to  have  higher  values  than  the  spectra  for  the  shorter  P-wave  windows. 
For  the  05/06/87  DM  event  the  difference  between  the  short-  and  long-window 
P-wave  spectra  is  seen  as  a  steady  decline  in  the  short-window  to  long-window 
ratio.  For  the  03/12/87  SR  event,  the  short-window  P  spectrum  is  most 
severely  depleted  relative  to  the  long-window  P  spectrum  at  specific  frequen¬ 
cies  (viz  between  1 .0  and  2.0  Hz). 

The  lower  right-hand  plots  in  the  figures  show  the  effects  on  the  Lg/P 
spectral  ratios.  As  would  be  expected  from  individual  spectral  observations, 
the  Lg/P  ratios  for  the  short  P-wave  windows  fall  below  those  obtained  for  the 
longer  P-wave  windows  over  most  of  the  frequency  band  from  1 .0  to  8.0  Hz. 
This  tendency  appears  to  be  the  same  for  events  at  both  the  SR  and  DM  test 
sites. 


32 


Spectral  analys 


Figure  13.  Spectral  analysis  of  the  WMQ  record  from  the  05/06/87  DM  explosion 


in  the  subsequent  analyses  we  have  used  only  the  longer  (viz  25 
seconds)  P-wave  windows.  The  main  effect  on  the  Lg/P  ratios  presented 
below  of  using  the  long  versus  the  short  window  is  expected  to  be  an  overall 
decrease  in  the  ratio.  Although  we  would  not  expect  this  effect  to  change  too 
drastically  from  event  to  event,  there  is  some  evidence  from  the  studies  just 
described  that  some  differences  with  frequency  could  occur  between  events  at 
different  sites.  This  problem  will  require  further  study. 

3.5  Lg/P  Spectral  Ratios 

Lg/P  spectral  ratios  were  computed  for  the  1 9  East  Kazakh  explosions  and 
27  regional  earthquakes  in  the  WMQ  sample.  Figures  14  through  18  show  the 
results  of  the  processing  performed  on  three  SR  explosions  and  two  regional 
earthquakes.  Comparing  the  spectra  for  the  explosions,  the  shapes  appear 
remarkably  consistent  over  the  magnitude  range  from  5.0  mb  to  6.1  mb  covered 
by  these  events.  In  all  cases  the  regional  P  spectra  show  a  relatively  flat  por¬ 
tion  extending  from  0.8  Hz  to  more  than  3  Hz.  Above  3  Hz  the  spectra  show  a 
steady,  rapid  decline.  The  Lg  spectra  for  the  same  explosions  show  a  less 
rapid  but  steady  decline  over  the  entire  frequency  range  plotted,  with  ampli¬ 
tudes  falling  below  the  regional  P  spectral  level  between  1.0  Hz  and  about  2.0 
Hz.  The  regional  P  and  Lg  spectra  for  the  two  earthquakes  also  have  roughly 
similar  shapes  (cf.  Figures  17  and  18)  even  though  the  two  events  are  at 
somewhat  different  ranges.  The  regional  P  spectra  again  show  a  nearly  flat 
portion  from  0.8  Hz  to  about  3.0  Hz  followed  by  a  rapid  decline  toward  higher 
frequencies.  The  Lg  spectra  show  a  steady  decline  starting  at  low  frequencies 
and  extending  to  4  Hz  or  so  before  decreasing  much  more  rapidly  at  higher 
frequencies.  In  contrast  to  the  explosions,  the  Lg  spectra  for  the  earthquakes 
do  not  reach  the  P-wave  spectral  level  until  a  frequency  of  5  Hz  or  more. 

Perhaps  more  interestingly,  comparing  the  earthquake  and  explosion  of 
equivalent  magnitude  (viz  the  12/14/86  and  the  06/14/88  events),  the  regional 
P-wave  spectra  for  the  two  events  match  very  closely  even  though  the  ranges 
are  again  quite  different.  In  contrast,  the  Lg  spectrum  for  the  earthquake  has 
much  more  high-frequency  energy.  It  needs  to  be  emphasized,  however,  that 
it  is  unclear  to  what  extent  attenuation  differences  because  of  the  different 
ranges  and  paths  may  affect  this  observation.  Next,  comparing  the  spectra  for 
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Shagan  Shot  at  WMQ,  Sept.  14,  1988,  04:01:37.59,  JVE 


Figure  14.  Spectral  analysis  of  the  WMQ  record  from  the  09/14/88  (JVE)  SR  oxplosi 


Shagon  Shot  ot  WMQ,  March  12,  1987,  01:57:17.2,  Mog.  5.5 
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the  5.3  mb  earthquake  (06/17/88)  and  the  5.5  mb  explosion  (03/12/87),  the 
regional  P-wave  spectra  are  again  seen  to  match  quite  closely  between  the 
two  events.  However,  the  earthquake  Lg  spectrum  is  roughly  a  factor  of  ten  or 
so  larger  than  the  corresponding  Lg  spectrum  from  the  explosion.  The  biggest 
difference  between  the  Lg  spectra  for  these  two  latter  events  appears  to  occur 
at  a  frequency  of  about  4  Hz.  In  general,  then,  the  Lg  spectra  for  the  regional 
earthquakes  appear  to  be  richer  in  high  frequencies  than  those  from  the  SR 
explosions  with  comparable  regional  P  spectra  recorded  at  WMQ. 

This  result  is  also  apparent  in  the  Lg/P  spectral  ratios  computed  for  these 
events  and  plotted  in  the  lower  right  portion  of  the  figures.  The  explosion 
spectral  ratios  are  observed  to  drop  off  precipitously  from  very  high  values  at 
frequencies  near  1  Hz  to  values  about  a  factor  of  40  lower  above  about  3  Hz. 
The  earthquake  Lg/P  spectral  ratios,  on  the  other  hand,  start  out  at  values  near 
ten  or  so  which  are  only  slightly  larger  than  the  explosion  ratios.  However,  the 
earthquake  ratios  decrease  by  only  about  a  factor  of  ten  within  the  frequency 
range  shown.  Thus,  these  results  indicate  that  Lg/P  spectral  ratios  for  the 
earthquakes  are  relatively  enriched  at  high  frequencies  compared  to  the  explo¬ 
sions  recorded  at  WMQ. 

The  Lg/P  spectral  ratios  from  the  1 1  Sh  explosions  were  averaged  and  the 
standard  deviation  computed  assuming  a  log-normal  distribution.  Similar 
averaging  was  performed  on  the  8  DM  explosions  and  the  27  regional  earth¬ 
quakes.  Figure  19  presents  the  Lg/P  spectral  ratio  averages  for  the  SR  explo¬ 
sions  and  the  earthquakes.  Differences  between  the  SR  and  DM  explosion 
ratios  will  be  described  in  the  following  section  of  this  report.  As  was  observed 
for  the  individual  ratios  presented  above,  the  average  Lg/P  ratios  for  earth¬ 
quakes  and  explosions  appear  to  be  indistinguishable  at  frequencies  near  1 
Hz.  However,  the  average  explosion  spectral  ratio  is  seen  to  drop  off  by  nearly 
two  orders  of  magnitude  toward  higher  frequencies  (  =  6-7  Hz).  Over  the  same 
interval  the  average  Lg/P  ratio  for  the  earthquakes  decreases  by  only  about 
one-half  order  of  magnitude.  As  a  result,  the  Lg/P  ratios  are  well  separated  at 
high  frequencies.  From  the  plots  of  the  mean  plus  and  minus  one  standard 
deviation  in  Figure  19,  we  see  complete  separation  of  the  earthquake  and 
explosion  Lg/P  spectral  ratios  above  about  2  Hz.  At  higher  frequencies  the 
ratios  are  separated  even  further. 
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These  observations  support  our  findings,  from  the  band-pass  measure¬ 
ments,  that  Lg/P  spectral  ratios  for  earthquakes  are  larger  at  high  frequencies 
than  for  similar  explosions  recorded  at  WMQ.  These  lates.  analyses  confirm 
the  same  results  for  larger  samples  of  both  explosions  and  earthquakes.  It 
should  also  be  noted  that  these  observations  appear  to  be  consistent  with  our 
previous  findings  (cf.  Murphy  and  Bennett,  1982;  Bennett  and  Murphy,  1986) 
for  nearby  explosions  and  earthquakes  in  the  western  U.3.  In  particular,  it 
was  found  there  that  the  earthquake  U  signals  were  relatively  enriched  at  high 
frequencies  while  the  regional  P  spectra  from  the  earthquakes  and  explosions 
appeared  to  have  quite  similar  shapes.  In  the  following  section  of  this  report, 
we  show  results  that  indicate  that  differences  in  the  Lg/P  ratios  between  SR 
and  DM  explosions  are  inadequate  to  diminish  the  apparent  discrimination 
capability  of  this  measurement.  The  one  outstanding  problem  appears  to  be 
the  effects  of  regional  path  differences  on  the  explosion  anH  earthquake  sig¬ 
nals  propagating  to  WMQ.  Since  we  cannot  identify  co-located  jarthquakes 
and  explosions  in  the  areas  of  interest,  the  best  approach  to  handling  these 
factors  appears  to  be  by  compensating  for  attenuation  differences.  To  achieve 
this  objective  a  better  understanding  of  attenuation  effects  on  regional  P  and  Lg 
signals  for  regions  of  interest  in  and  around  the  Soviet  Union  will  be  required. 
Results  of  some  initial  investigations  toward  this  objective  are  presented  in 
Section  V  of  this  report.  In  future  work  we  will  be  performing  similar  analyses 
on  regional  signals  recorded  at  other  high  quality  digital  stations,  including  the 
Soviet  IRIS  stations,  and  trying  to  develop  a  better  theoretical  understanding  of 
the  Lg/P  spectral  ratio  behavior. 
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IV.  Comparison  of  Shagan  River  and  Degelen 
Mountain  Explosions  Recorded  at  WMQ 


4.1  SR  and  DM  Bias  Differences 

Application  of  regional  phase  signal  measurements  to  discrimination  and 
yield  estimation  has  been  hampered  over  the  years  by  the  lack  of  theoretical 
understanding  of  he  excitation  of  these  signals  by  explosion  and  earthquake 
sources.  Of  particular  interest  in  this  regard  is  whether  different  test  sites  pro¬ 
duce  different  regional  phase  signals  from  otherwise  comparable  events.  In 
proposing  that  Lg  magnitudes  could  be  used  to  estimate  the  yields  of  under¬ 
ground  nuclear  explosions,  Nuttli  (1986a,b)  argued  that  explosions  with 
equivalent  energy  release  (i.e.  yield)  should  produce  Lg  signals  with  equal 
amplitudes  near  the  source.  Thus,  the  Lg  magnitude  should  provide  a  measure 
of  explosion  yield  independent  of  the  test  site  where  the  event  occurred.  Nuttli 
based  this  conclusion  cn  a  small  number  of  explosions  with  announced  yields 
in  hard  rock  at  a  few  diffe.ent  test  sites. 

Nuttli  further  suggested  that  observations  of  the  differences  between  Lg 
magnitude  and  teleseismic  mb  from  different  test  sites  could  be  related  to 
attenuation  of  the  teleseismic  P  waves  in  the  upper  mantle  under  the  test  sites. 
Thus,  the  difference  between  Lg  magnitude  and  mb  should  measure  the  telese¬ 
ismic  mb  bias  of  the  test  site.  Nuttli  (1986a,b)  determined  the  test  site  biases 
for  NTS  and  SR  to  be  -0.31  and  0.036,  respectively.  Furthermore,  Nuttli 
(1987)  estimated  the  bias  for  DM  to  be  0.27.  The  large  difference  in  biases 
between  SR  and  DM  was  surprising  considering  that  the  two  test  sites  are 
separated  by  only  about  50  km.  Based  on  this  result  Nuttli  (1987)  concluded 
that  either  P-wave  amplitudes  from  DM  explosions  are  1.9  times  larger  than 
SR  explosions  or  Lg  amplitudes  are  1.9  times  smaller  for  DM  explosions.  To 
obtain  additional  insight  into  this  problem,  we  have  analyzed  the  behavior  of 
the  regional  P  and  Lg  signals  recorded  at  WMQ  from  SR  and  DM  events. 
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4.2  Regional  Waveform  Comparisons 


As  noted  above  in  Section  II,  our  database  of  explosions  recorded  at 
WMQ  includes  18  SR  explosions  and  9  DM  explosions.  We  presented  the 
vertical  component  records  at  WMQ  for  all  of  these  events  above  in  Figure  2. 
Figure  20  shows  three  examples  from  each  test  site  on  a  somewhat  expanded 
scale  for  comparison.  The  selected  events  from  each  test  site  have  roughly 
comparable  magnitudes  between  4.8  and  5.9  mb.  The  records  in  Figure  20 
appear  to  show  nearly  the  opposite  behavior  to  that  observed  by  Nuttli.  In  par¬ 
ticular,  the  regional  P  phases  from  the  DM  explosions  are  in  all  cases  smaller 
than  the  Lg,  whereas  the  regional  P  phases  from  the  SR  explosions  are  larger 
than  the  Lg.  One  of  the  main  differences  appears  to  be  the  large,  relatively 
high-frequency,  initial  P  phase  for  the  SR  explosions  which  is  not  apparent  for 
the  DM  explosions.  One  explanation,  consistent  with  Nuttli’s  observations, 
might  be  that  for  SR  explosions  some  of  the  high-frequency  energy  normally 
contributing  to  teleseismic  P  is  being  trapped  in  the  regional  P  waveguide;  but 
several  other  possible  explanations  also  need  to  be  evaluated. 

Figure  21  shows  comparisons  of  the  regional  P  and  Lg  signal  levels  at  the 
two  East  Kazakh  test  sites  for  all  1 8  SR  explosions  and  the  9  DM  explosions 
as  measured  at  WMQ.  The  plot  presents  the  logarithms  of  the  maximum  P 
amplitudes  measured  from  the  broadband  records  versus  the  logarithms  of  the 
maximum  Lg  amplitudes  also  measured  from  the  broadband  records.  The  data 
show  considerable  scatter  and  some  intermingling  but  are  separated  on  aver¬ 
age.  Regression  lines  were  fit  to  each  data  set.  These  are  described  by 

log  P(max)  =  0.744  log  Lg(max)  +  0.909 


for  SR  and 


log  P(max)  =  0.899  log  Lg(max)  +  0.147 

for  DM.  Within  the  range  of  overlap  of  the  Lg  amplitude  measurements,  the 
maximum  regional  P-wave  amplitudes  are  separated  by  from  0.2  to  0.4 
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Figure  21. 


Regional  P  versus  Lg  maximum  amplitude  levels 
for  SR  and  DM  explosions  recorded  at  WMQ. 
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magnitude  units  on  average  between  the  two  test  sites.  As  with  the  examples 
shown  above,  the  regional  P  phases  from  the  SR  explosions  for  the  total  sam¬ 
ple  appear  to  be  on  average  a  factor  of  about  two  larger  in  amplitude  than 
regional  P  phases  from  DM  explosions  which  have  comparable  Lg  amplitudes 
at  WMQ.  It  should  be  noted  that  the  standard  deviations  of  the  measurements 
relative  to  these  regressions  are  large.  It  is  also  interesting  that  the  regression 
lines  tend  to  converge  at  larger  amplitudes  (i.e.  larger  magnitude),  but  this  is  to 
a  large  measure  controlled  by  the  single  larger  amplitude  observation  at  DM 
which  is  intermingled  with  the  SR  data  points. 

4.3  Spectral  Comparisons 

We  computed  spectra  for  the  initial  P  and  Lg  windows  at  WMQ  from  com¬ 
parable  magnitude  SR  and  DM  explosions,  whose  time  histories  were 
presented  above  in  Figure  20.  These  are  shown  in  Figure  22.  In  computing 
these  spectra  we  used  a  short  window  (viz  12.8  seconds)  including  the  initial  P 
and  a  long  window  (viz  102.4  seconds)  which  encompassed  nearly  all  the  Lg 
energy.  The  time  histories  were  tapered  and  the  spectra  smoothed  using  a 
running  average.  The  instrument  response  for  the  WMQ  broadband  instrument 
was  also  removed  prior  to  plotting  the  spectra  in  Figure  22.  Aside  from  some 
minor  fluctuations  at  certain  frequencies,  the  Lg  spectra  appear  to  be  quite  con¬ 
sistent  between  the  comparable  SR  and  DM  •'xplosions.  Although  the  Lg  spec¬ 
tra  for  the  DM  explosions  in  both  cases  tern  slightly  below  those  for  the  SR 
explosions  at  a  frequency  near  1  Hz,  this  difference  might  be  explained  by  the 
slightly  lower  magnitudes  of  the  DM  events.  The  differences  in  the  regional  P- 
wave  spectra  between  the  SR  and  DM  events  cover  a  broader  range  of  fre¬ 
quencies  and  are  generally  larger.  The  latter  differences  probably  cannot  be 
explained  by  the  small  magnitude  differences  between  the  compared  events. 
The  observations  suggest  that  the  regional  P  waves  from  the  SR  explosions 
are  enhanced  over  a  broad  range  of  frequencies  relative  to  DM  explosions  with 
comparable  Lg  signals. 

Spectral  analyses  were  also  performed  on  a  larger  sample  of  1 1  SR  and  8 
DM  explosions  using  the  same  procedures  as  those  described  in  Section  III  for 
determining  Lg/P  spectral  ratios.  Figures  23  through  26  show  the  results  of 
these  analyses  performed  on  two  SR  and  two  DM  underground  nuclear 
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Figure  22.  Spectral  comparisons  for  initial  P  and  Lg  phases  from 
similar  magnitude  explosions  at  SR  and  DM  recorded  at 
WMQ.  Note  that  instrument  response  has  been  removed 
from  spectral  shape. 
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Figure  23.  Spectral  analysis  of  the  WMQ  record  from  the  06/14/88  SR  explosion. 


06.7,  Mog.  5.3 


Figure  2C> .  Spectral  analysis  of  the  WMQ  record  from  the  1  1/2  3/88  DM  explosion. 


explosions  with  similar  magnitudes  in  the  range  5.0  to  5.5  mb.  The  ground 
motion  records,  shown  at  the  top  of  the  figures,  indicate  the  same  tendency 
noted  above  with  relatively  larger  regional  P  phases  for  the  SR  events  versus 
DM  events  with  comparable  Lg  signals.  The  results  presented  for  the  regional 
P  spectra  in  Figures  23  through  26  correspond  to  a  longer  P-wave  window  (viz 
25  seconds  long  and  beginning  1  second  prior  to  the  initial  P  markers  shown 
superimposed  on  the  waveform  traces).  The  Lg  windows  covered  a  group 
velocity  range  from  3.5  to  2.9  km/sec,  indicated  by  the  later  markers,  and 
clearly  include  most  of  the  Lg  energy.  Although  some  additional  Lg  energy  is 
seen  to  arrive  after  the  end  of  the  defined  windows,  we  experimented  with  win¬ 
dows  which  included  this  energy  also  and  found  that  the  choice  of  the  Lg  win¬ 
dow  bounds  made  little  difference  in  the  spectral  shape  or  level.  The  windows 
were  tapered;  discrete  Fourier  transforms  were  then  computed  and  resampled 
and  smoothed  to  provide  equivalent  frequency  spacing  for  the  regional  P  and 
Lg  signal  spectra.  Finally,  Lg/P  spectral  ratios  were  computed  at  each  fre¬ 
quency. 

The  spectra  in  Figures  23  to  26  appear  remarkably  similar  from  event  to 
event.  In  particular,  over  the  frequency  band  shown  in  these  figures,  the  Lg 
spectral  shapes  can  be  made  to  overlay  almost  exactly  between  all  four 
events.  Slight  differences  in  the  general  level  of  the  Lg  spectra,  however,  are 
apparent  from  event  to  event.  The  regional  P  spectra  are  slightly  more  vari¬ 
able  in  both  shape  and  level.  Comparing  first  the  two  DM  shots,  which  are 
equal  in  magnitude  (viz  5.3  mb),  the  only  real  difference  is  a  very  slight 
increase  in  the  regional  P-wave  spectral  level  for  the  06/06/87  event  over  the 
11/23/88  event  in  the  limited  frequency  band  from  about  2.0  to  5.0  Hz  (cf.  Fig¬ 
ures  23  and  24).  Comparing  next  the  06/06/87  DM  event  with  the  03/12/87  SR 
event  (5.5  mb),  we  see  that  the  Lg  spectra  match  very  closely  in  both  shape 
and  level  (cf.  Figures  23  and  25)  even  though  the  teleseismic  mb  magnitudes 
are  different  by  0.2  magnitude  units.  In  contrast,  the  regional  P-wave  spectrum 
for  the  03/12/87  SR  event  lies  generally  above  that  of  the  06/06/87  DM  event 
with  the  greatest  differences  (up  to  a  factor  of  three  in  amplitude)  occurring  in 
the  frequency  range  from  1 .0  to  3.0  Hz.  Finally,  comparison  of  the  06/06/87 
DM  event  with  the  06/14/87  SR  event  (5.0  mb)  indicates  that  the  regional  P- 
wave  spectra  match  more  closely  in  shape  and  level  (cf.  Figures  23  and  26); 
while  the  Lg  spectrum  for  the  smaller  SR  event  lies  somewhat  below  that  of  the 
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DM  event  over  a  limited  frequency  band  from  about  2.0  to  5.0  Hz.  These 
resuits  again  appear  to  be  consistent  with  relatively  greater  excitation  of 
regional  P  signals  at  some  frequencies  from  SR  explosions  which  have  the 
same  level  of  Lg  excitation  as  DM  explosions. 

The  latter  observation  shows  up  as  a  relatively  subtle  effect  in  the  Lg/P 
spectral  ratios  presented  in  Figures  23  to  26.  Comparing  first  the  Lg/P  ratios 
for  the  shots  at  common  test  sites,  it  is  seen  that  the  ratios  match  quite  closely 
over  the  entire  frequency  band  presented.  Comparisons  of  the  LJP  spectral 
ratios  between  SR  and  DM  events  indicates  that  the  SR  spectra  generally  lie 
slightly  below  the  DM  spectra  at  almost  all  frequencies.  Figure  27  shows  this 
same  comparison  for  the  suite  of  1 1  SR  events  and  8  DM  events.  The  figure 
presents  the  average  Lg/P  spectral  ratios  from  these  events  along  with  ±  1(o) 
bounds  on  these  mean  values  as  a  function  of  frequency.  The  L^/P  ratios  are 
observed  to  be  nearly  identical  around  1  Hz  and  more  separated  at  higher  fre¬ 
quencies  with  the  SR  average  lying  below  the  DM  average.  Over  most  of  the 
frequency  band  from  1  Hz  to  8  Hz,  the  separation  appears  to  hold  at  the  ±  1  (o) 
level,  although  just  barely. 

From  these  analyses  we  would  conclude  that  SR  sources  tend  to  be  rela¬ 
tively  more  efficient  in  generating  regional  P  than  DM  sources  with  equivalent 
Lg  signals.  Above  we  suggested  an  explanation  consistent  with  Nuttli’s  obser¬ 
vation  of  larger  teleseismic  P  from  DM  events  might  be  preferential  trapping  of 
P  energy  from  the  SR  explosions  in  the  crust  and  upper-mantle  waveguide. 
Alternatively,  these  observations  of  relatively  larger  Lg/P  ratios  from  DM  events 
could  also  be  interpreted  as  greater  excitation  of  Lg  by  the  DM  sources.  The 
latter  might  be  caused  by  more  efficient  coupling  between  the  Lg  crustal 
waveguide  at  the  DM  source  (in  contrast  to  SR)  and  the  transmission  path 
waveguide  to  WMQ.  If  such  effects  exist  they  have  not  been  demonstrated, 
and  additional  explanation  of  Nuttli’s  teleseismic  P  differences  would  be 
required.  Other  explanations  for  these  observations  also  exist  and  require 
further  study.  We  plan  follow-on  studies  using  the  Soviet  IRIS  and  other 
regional  stations  to  further  resolve  these  differences  between  SR  and  DM. 
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It  should  also  be  noted  that  the  differences  between  SR  and  DM  tend  to 
make  the  Lg/P  spectral  ratios  for  the  DM  events  appear  slightly  more 
earthquake-like  (cf.  Section  III  above).  However,  it  can  be  seen  by  comparing 
the  observations  in  Figure  27  with  the  discrimination  results,  presented  above 
in  Figure  19,  that  there  still  appears  to  be  good  separation  between  the  earth¬ 
quake  and  explosion  sample  above  about  3  Hz.  We,  therefore,  conclude  that 
test  site  differences  like  those  between  SR  and  DM  are  not  sufficient  to  con¬ 
found  the  discrimination  differences  seen  in  the  Lg/P  spectral  ratios  between 
underground  nuclear  explosions  and  earthquakes. 


V.  Investigation  of  Path  Effects  on  Regional  Signals 


5.1  Attenuation  Measurements 

As  noted  above,  one  argument  against  the  effectiveness  of  empirically 
derived  discriminants,  such  as  the  Lq/P  spectral  ratios  presented  in  Section  III, 
is  that  they  do  not  take  into  account  the  differences  in  propagation  path  effects 
when  sources  are  not  co  located.  The  results  we  showed  above  for  East 
Kazakh  explosions  compared  the  observed  ratios  to  those  from  earthquakes, 
which  in  some  cases  had  similar  ranges  to  WMQ,  but  whose  epicenters  were 
located  in  many  cases  along  the  southern  Soviet  border  (cf  Figure  3  above) 
t  herefore,  the  propagation  paths  to  the  station  from  underground  nuclear  tests 
are  frequently  quite  different  than  those  from  the  earthquakes  used  for  com 
pansons.  In  our  analysis  we  attempted  to  compensate  for  these  effects  to 
some  degree  by  including  earthquakes  from  a  variety  of  locations  so  as  not  to 
bias  the  result  toward  a  particular  propagation  path  from  a  single  earthquake 
source  area  However,  a  clearly  betler  approach  would  be  to  handle  the  pro 
pagatiori  effects  more  directly  by  making  appropriate  adjustments  to  the  signal 
or  signal  spectra  from  each  source  area  for  attenuation  differences  Develop 
merit  of  a  better  understanding  of  regional  propagation  effects  for  paths  of 
interest  in  the  Soviet  Union  is  a  major  task  We  have  attempted  to  make  some 
progress  on  this  problem  during  the  past  year  by  analyzing  I  attenuation  for 
paths  to  selected  GD.'.N  and  Soviet  IRIS  stations  from  the  East  Kazakh  source 
region  In  the  future  it  is  anticipated  that  similar  information  can  be  developed 
into  a  regionalized  attenuation  model  for  the  Soviet  Union  and  surrounding 
areas  of  F  urasia  which  can  in  turn  tie  used  to  correct  regional  phase  spectra 
for  propagation  effects  and  identify  source  related  differences 

l!  is  clear  that  tegional  phase  signals  are  affected  by  attenuation  In  addi 
turn  to  the  normal  amplitude  decrease  associated  with  geometrical  spreading. 
lh’  regional  phase  signal  loses  energy  due  to  inherent  medium  absorption  a  rid 
scattering  try  inhomogeneities  m  the  waveguide  In  developing  the  |  magni 
hide  scale.  Nutlli  (I'ict)  used  an  analytic  expression  tor  geometrical  spreading 
and  an  effec  live  attenuation  teim.  including  both  absorption  and  scattering 


ai 


effects,  to  determine  the  near-source  amplitude  level.  According  to  this  model 
of  the  Lg  signal,  propagating  energy  is  dissipated  in  relation  to  the  time  it 
spends  in  the  medium,  with  the  higher  frequencies  being  preferentially 
absorbed.  Since  tha  time  a  number  of  methods  have  evolved  for  deriving  the 
effective  Q,  representative  of  attenuation,  from  observational  data.  Some  of 
these  methods  involve  measurement  of  the  coda  decay  (e.g.  Aki  and  Chouet, 
1975;  Herrmann,  1980;  Nuttli,  1981,  1986a,b)  while  others  are  based  on  obser¬ 
vation  of  the  behavior  of  the  signal  spectrum  at  one  or  more  stations  (e.g. 
Chun  et  at.,  1987;  Baumgardt  and  Ziegler,  1988;  Sereno  et  at.,  1988).  Infor¬ 
mation  on  coda  behavior,  spectral  content,  and  calibration  events  can  be  used 
to  derive  an  effective  Q0  for  the  path  from  SR  to  WMQ.  In  the  current  study 
we  have  used  the  Lg  amplitude  differences  between  WMQ  and  HIA  and  the 
Soviet  IRIS  stations  for  common  SR  explosions  to  derive  Q0’s  for  those  paths 
as  well. 

5.2  Lg  Attenuation  for  the  Path  to  WMQ 

In  our  initial  analysis  of  Lg  attenuation,  we  have  focused  on  the  path  to 
station  WMQ  from  SR  events.  As  noted  above,  a  relatively  large  database  of 
strong  Lg  signals  exists  for  this  path.  In  a  previous  report  Bennett  et  at.  (1989) 
noted  that  the  RSTN  station  RSSD  recorded  strong  Lg  signals  from  NTS  explo¬ 
sions  and  that  the  general  characteristics  of  the  regional  signals  at  WMQ  and 
RSSD  appeared  quite  similar,  even  though  the  stations  were  located  at  some¬ 
what  different  distances  from  the  respective  source  areas,  960  km  from  SR  to 
WMQ  and  1270  km  from  NTS  to  RSSD.  To  develop  the  effective  attenuation 
for  the  paths  to  these  stations,  we  have  compared  the  amplitude  and  spectral 
characteristics  of  the  Lg  signals  from  explosions  of  equivalent  yield  at  NTS  and 
SR  to  determine  to  what  extent  those  characteristics  can  be  attributed  to 
attenuation  differences. 

We  first  selected  several,  large  magnitude  events  which  were  well- 
recorded  at  each  station.  Vertical  component  records  for  these  events  are 
shown  in  Figure  28.  They  include  three  SR  explosions  recorded  at  WMQ  and 
four  NTS  explosions  recorded  at  RSSD.  Table  5  summarizes  the  source  infor¬ 
mation  on  the  seven  events.  To  make  the  signal  records  more  comparable 
between  the  two  stations,  we  removed  the  different  instrument  responses  from 
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Table  5.  Source  Information  for  Similar 
Yield  Explosions  at  NTS  and  SR 


Event 

Lat. 

Lon. 

mb 

NTS 

03/01/84 

37.07 

-116.05 

5.9 

05/31/84 

37.10 

-116.05 

5.8 

05/02/85 

37.25 

-116.33 

5.7 

12/05/85 

37.05 

-116.05 

5.7 

SR 

12/13/87 

49.99 

78.84 

6.1 

04/03/88 

49.92 

78.95 

6.1 

09/14/88 

49.83 

78.81 

6.1 
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the  original  records  and  reconvolved  with  the  instrument  response  for  the 
short-period  WWSSN  station  to  produce  the  records  in  Figure  28.  The  approx¬ 
imate  beginning  of  the  Lg  windows  is  designated  by  the  arrows  at  a  group  velo¬ 
city  of  3.6  km/sec. 

As  noted  above  in  Section  IV,  Nuttli  (I986a,b)  estimated  the  teleseismic 
mb  bias  for  NTS  to  be  -0.31  and  for  SR  to  be  0.036.  Thus,  equivalent-yield 
explosions  at  NTS  and  SR  would  be  expected  to  have  mb  values  about  0.35 
magnitude  units  larger  for  the  SR  explosions.  In  selecting  the  explosions  in 
Table  5  we  attempted  to  pick  events  with  approximately  this  magnitude 
difference,  so  that  the  explosions  at  the  two  test  sites  are  nearly  equivalent  in 
yield.  In  particular,  pairs  of  explosions  of  nearly  equivalent  yield  from  each  test 
site  were  identified  and  Lg  spectra  determined  from  the  time  histories  recorded 
at  the  respective  stations.  Figure  29  shows  the  Lg  spectra  for  two  such  explo¬ 
sion  pairs  (viz  the  04/03/88  SR  event  paired  with  the  05/31/84  NTS  event,  on 
the  left,  and  the  09/14/88  SR  event  paired  with  the  05/02/85  NTS  event,  on  the 
right).  The  spectra  plotted  in  the  figure  include  the  common  WWSSN  response 
referred  to  above.  The  Lg  spectra  are  observed  to  be  fairly  sharply  peaked  at 
a  frequency  between  0.5  Hz  and  1.0  Hz.  The  sharpness  of  the  peak  is  to 
some  extent  related  to  the  WWSSN  response;  the  broadband  spectra 
described  above  in  Sections  II  and  111  had  less  exaggerated  peaks.  The  spec¬ 
tral  shapes  appear  to  be  consistent  with  the  expected  behavior  of  increased 
attenuation  causing  depletion  of  high-frequency  energy  at  the  more-distant  sta¬ 
tion.  The  Lg  spectra  at  WMQ  show  a  steady  decline  above  the  peak  before 
reaching  an  apparent  noise  floor  at  about  6  Hz.  The  Lg  spectra  at  RSSD 
decline  more  rapidly  above  the  peak  and  reach  the  apparent  noise  level  at  a 
much  lower  frequency  (viz  about  3  Hz).  Lg  spectral  ratios  determined  for  the 
two  event  pairs  over  a  frequency  band  from  0  to  3  Hz  are  presented  in  Figure 
30.  The  measured  Lg  spectral  ratios  are  remarkably  consistent  between  the 
two  pairs.  Ignoring  the  variations  at  very  low  frequencies,  which  may  be 
related  to  inaccurate  response  information,  the  Lg  spectral  ratio  (WMQ  to 
RSSD)  shows  a  steady  increase  over  the  range  from  0.2  to  3.0  Hz.  A  slight 
tendency  to  roll-over  and  flatten-out  toward  the  higher-frequency  end  is  also 
apparent. 
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WMQ/RSSO  SPECTRA 


Figure  30.  Lg  spectral  ratios  WMQ/RSSD  for  event  pairs 
from  Figure  29. 
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In  analyzing  this  observation  it  is  assumed  that  the  observed  Lg  spectrum 
from  an  explosion  source  can  be  represented  as 

A(f)  =  A0(f)  1(f)  G(R)  S(f)  H(f,R)  (1) 

where  A0  is  the  source  function  effect,  I  is  the  instrument  response,  G  is 
geometrical  spreading,  S  is  site  response  and  H  is  path  attenuation.  Assuming 
equivalent  source  functions,  instrument  response,  and  station  effects,  then  the 
spectral  ratio  of  the  signals  for  the  event  pairs  can  be  written  as 

A,(f)  G(R,)  H, (f,  Ri) 

P(f)  =  — H -  =  — 1— 11  —Hi — H  (2) 

A2(f)  G(R2)  H2(f,  R2)  1  ' 

where 

m/i  rcfAt  Jif  R  /o\ 

H(f.R)  =  exp  — ~  =  exp  — ■  (3) 

and 


Thus, 


So,  the  path  Q  functions  can  be  used  to  predict  the  spectral  ratios  for  the  event 
pairs. 


In  these  analyses  we  took  the  common  representation  of  O  as 


0  =  00^ 


(6) 


Studies  of  Lg  attenuation  in  continental  areas  have  typically  found  the  fre¬ 
quency  exponent,  c,  to  have  values  between  0.4  and  0.6  (e.g.  Mitchell,  1990; 
Sereno  et  al.,  1988)  In  the  present  study  we  have  not  attempted  to  determine 
the  frequency  dependence  of  Q  but  havp  instead  assumed  values  in  the 
appropriate  range.  Therefore,  we  initially  assumed,  based  on  preliminary 
analysis,  that  c  was  0.4  for  the  path  from  SR  to  V.MQ  and  0.5  for  the  path 
from  NTS  to  RSSD.  Using  this  as  a  constraint,  we  determined  values  of  Q0 
which  would  give  the  appropriate  spectral  ratios  when  applied  in  Equation  5. 
Figure  31  shows  the  predicted  spectral  ratios  for  the  explosion  pairs  under  the 
assumption  that  Q  =  285  f0  4  for  the  path  to  WMQ  and  Q  =  200  f°  5  for  the  path 
to  RSSD.  The  agreement  between  the  predicted  and  observed  behavior  for 
the  Lg  spectral  ratio  appears  to  be  reasonable,  but  the  estimates  of  Q  could  be 
refined  to  provide  even  closer  agreement.  Furthermore,  it  should  be  noted 
that,  even  though  the  observed  Lg  spectral  ratios  are  reasonably  matched  by 
predictions  based  on  the  path  Q’s,  these  Q  values  do  not  produce  Lg  magni¬ 
tudes  which  would  give  the  right  teleseismic  mb  bias  with  respect  to  mb(Lg).  In 
particular,  these  Q  values  give  Lg  magnitudes  which  are  too  large  for  both 
paths. 

The  magnitude  bias  and  the  spectral  ratio  results  can  be  made  consistent 
by  adjusting  the  attenuation  for  both  paths.  We,  therefore,  proceeded  under 
the  constraint  that  the  Lg  magnitudes  measured  from  the  records  in  Figure  28 
should  give  the  appropriate  test  site  bias,  as  determined  by  Nuttli  (1986a,b), 
when  compared  to  the  reported  teleseismic  mb.  The  logarithms  of  the  peak  Lg 
amplitudes  measured  from  the  records  in  Figure  28  were  used  in  an  inversion 
scheme  to  determine  the  Q0  values  for  each  event  which  would  give  Lg  magni¬ 
tudes  that  provided  the  correct  test  site  bias  (viz  mb  -  m0(Lg)  =  -0.31  for  NTS 
and  mb  -  mb(Lg)  =  0.036  for  SR).  For  the  three  SR  events  recorded  at  WMQ  in 
Table  5,  we  determined  Q0  to  be  452  with  o  =  34.  The  four  NTS  events 
recorded  at  RSSD  (cf.  Table  5)  gave  a  Q0  of  239  with  a  =  5.  The  relatively 
small  o  values  suggest  that  the  bias  differences  were  quite  consistent  between 
events  recorded  at  these  stations. 
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I-’L<jure  31.  Comparison  of  observed  and  predicted  I,(J  sj 
tral  ratios  derived  from  attenuation  model 


As  a  final  check  on  these  effective  Q  values  for  the  paths  SR/WMQ  and 
NTS/RSSD,  we  reassessed  the  Lg  spectrai  ratio  for  the  comparable-yield 
explosions.  For  the  SR/WMQ  path  we  used  Q  =  452  f0-5,  and  for  the 
NTS/RSSD  path  we  used  Q  =  239  f0  5.  Equation  5  above  was  again  used  to 
predict  the  spectral  ratios.  It  can  be  seen  in  Figure  32,  that  these  attenuation 
models  do  an  excellent  job  of  predicting  the  observed  ratios.  It’s  possible  that 
attenuation  models  with  different  frequency  dependence  might  also  do  a  rea¬ 
sonable  job  of  fitting  these  observations.  However,  we  would  expect  that  only 
minor  variations  for  the  Q0’s  can  be  permitted  to  simultaneously  satisfy  the  test 
site  bias  constraints.  A  more  thorough  analysis  of  the  sensitivity  of  the 
attenuation  estimates  to  these  measurements  would  be  useful. 

5.3  Attenuation  for  Paths  to  Other  Regional  Stations 

As  pointed  out  above  in  Section  II,  strong  Lg  signals  have  been  recorded 
at  a  regional  network  of  stations  from  East  Kazakh  explosions.  In  this  section 
we  will  determine  the  relative  amplitudes  of  these  Lg  signals  between  stations 
and  use  ‘hese  observations  to  estimate  effective  Q0  values  for  the  paths  from 
SR  to  the  various  stations.  We  start  with  Nuttli’s  definition  of  the  Lg  magnitude: 


mb(Lg)  =  log  A(10km)  +  2.96 


where 


(7) 


A(10km)  =  A(R)  (^-)3 


sin(R/1 11.1) 
sin(  1 0/1 11.1) 


9 


exp 


Ttf(R-IQ) 
Q  U 


(8) 


where  A(R)  is  the  Lg  ground  motion  amplitude  in  pm  measured  on  a  WWSSN 
short-period,  vertical-component  instrument.  The  Lg  magnitude  difference 
between  two  stations  recording  the  same  event  should,  in  the  absence  of 
measurements  errors,  be  zero.  So,  we  can  write 


=  Amb'(Lg)  +  log 


Amb(Lg) 


exp 


ref 


=  0 

/ Ri  - 
\  U1Q1 


R2  -  10 

u2  Q2 


(9) 
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where 


A(R1)  -j  R,  i  sin(R1/1 11.1) 

-  109  air?  +  3  109  iq  +  I  109  ii^HTTTT  1101 


Then 


Ttf 


/  Ro  -  1 0  Ri  -  10  \ 

("uToT  -  1096  -  Am»  :L^) 


(11) 


and 


Qp  - 


T2  Q-i  n  f  log  e 


Ti  n  f  log  e  +  Qi  Amb'(Lg) 


(12) 


where  T  =  R/U  and  R»10.  In  these  analyses  the  U  amplitude  measurements 
were  made  from  a  seismogram  which  had  been  converted  to  WWSSN 
response,  so  we  assumed  that  f  =  1  Hz  and  that  the  and  Q2  values  then 
correspond  to  Q0's  for  the  station  involved.  Since  we  had  some  prior 
knowledge  of  the  attenuation  to  station  WMQ  (based  on  the  analyses 
presented  in  the  preceding  section)  and  a  iarge  amount  of  data  recorded  there, 
it  was  used  as  the  reference  station.  Lg  magnitude  residuals  relative  to  WMQ 
were  then  analyzed  for  the  Soviet  IRIS  stations  (viz  GAR,  ARU,  OBN,  and  KIV) 
and  for  the  other  useful  CDSN  station  (viz  HIA)  to  determine  the  respective  Q0 
values  for  these  paths  from  SR  events. 


The  first  step  in  this  process  was  t^  convert  all  instrument  responses  to  an 
equivalent  WWSSN  short-period  instrument.  To  accomplish  this  the  individual 
station  responses  were  deconvolved  from  the  observed  records  and  the  output 
reconvolved  with  a  standard  WWSSN  response.  Figures  33  and  34  illustrate 
the  application  of  this  process  to  the  WMQ  and  ARU  records  for  two  SR  explo¬ 
sions.  For  both  stations  the  conversion  process  is  seen  to  basically  remove 
some  of  the  lower-frequency  energy  in  the  records.  In  particular,  this  greatly 
reduces  the  low-frequency  microseismic  noise  (enhancing  the  apparent  signal- 
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figure  33 


to-noise)  and  also  reduces  the  low-frequency  surface  waves. 


We  next  measured  the  sustained  Lg  amplitudes  using  the  procedures 
described  by  Nuttli  (1973;  1986a,b)  and  computed  the  residuals  (i.e.  Amb'(Lg) 
for  the  common  events.  Table  6  summarizes  the  data  that  were  analyzed  in 
this  way.  The  residuals  for  each  event/station  path  were  converted  to  esti¬ 
mates  of  Q0  for  these  paths  using  the  relationship  in  equation  12.  The  average 
Q0  values  for  the  paths  are  shown  in  Figure  35.  The  values  range  from  a  low- 
Q  path  value  of  428  to  station  GAR  up  to  a  high-G  path  of  761  to  OBN.  The 
Q0  value  of  580  determined  for  the  path  to  KIV  is  surprisingly  large  considering 
the  low  Lg  signal  level  normally  observed  there  and  probably  represents  a  max¬ 
imum  estimate.  For  comparison  we  show  in  Figure  36  Q0  values  derived  by 
Nuttli  for  paths  to  Eurasian  WWSSN  stations  from  East  Kazakh  test  site 
events.  The  Q0  values  determined  in  the  present  study  are  quite  consistent 
with  Nuttli’s  values  for  similar  paths.  They  also  agree  generally  with  Q  values 
recently  reported  by  Xie  and  Mitchell  (1990)  based  on  spectral  studies  of  Lg  for 
central  Asian  paths  from  explosions  and  earthquakes.  In  particular,  we  find 
relatively  low  Q’s  to  the  south  of  the  test  site  consistent  with  high  attenuation  of 
Lg  in  the  complex  tectonic  region  of  the  southern  Soviet  border.  On  the  other 
hand,  Q’s  for  the  northwestern  paths  are  high  corresponding  to  low  attenuation 
of  Lg  in  the  stable  continental  platform  area  of  the  Soviet  interior.  The  largest 
uncertainties  on  the  Q0  values  estimated  from  the  Lg  magnitude  residuals  are 
associated  with  station  ARU  and  OBN.  We  would  expect  some  refinement  of 
these  estimates  as  additional  data  become  available.  We  also  plan  to  more 
fully  investigate  the  frequency  dependence  of  Q  over  the  rather  limited  band  of 
frequencies  which  these  data  will  permit. 


Table  6.  Q0  Estimates  Derived  From  Measured 
Lg  Magnitude  Residuals 


Date 

WMQ 

HIA 

ARU 

GAR 

OBN 

KIV 

09/14/88 

452 

559 

707 

433 

818 

571 

1 0/1 8/88 

452 

- 

641 

410 

- 

- 

11/12/88 

452 

- 

675 

451 

- 

- 

1 2/1 7/88 

452 

577 

645 

- 

734 

- 

02/1 2/89 

452 

560 

577 

- 

- 

- 

02/17/89 

452 

- 

507 

417 

- 

- 

07/08/89 

452 

575 

604 

428 

731 

589 

452* 

568±1 0 

622±66 

428±1 6 

761 ±49 

580±1 3 

*  The  Q0  value  at  the  reference  station  WMQ  was  fixed  at  452 
based  on  previous  analyses. 


Figure  35 


Path  Q0  values  to  IRIS  and  CDSN  stations  derived  from 
Lg  magnitude  residuals. 


ure  36.  Q0  values  determined  by  Nuttli  (1986b)  for  paths 
from  SR  to  WWSSN  stations. 
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VI.  Summary  and  Conclusions 


6.1  Summary  of  Research  Program 

The  research  described  here  represents  the  first  year’s  effort  of  a  three 
year  program  designed  to  evaluate  the  capability  of  regional  seismic  stations  to 
discriminate  underground  nuclear  explosions  from  earthquakes  and  non¬ 
nuclear  explosions.  This  effort  to  date  has  focused  on  investigation  of  the 
regional  signals  recorded  at  CDSN  and  Soviet  IRIS  stations.  These  two 
sources  of  data  have  become  available  only  within  the  last  few  years  and  an 
experience  base  needs  to  be  developed  in  order  to  understand  their  capability 
for  monitoring  seismic  events  within  the  Soviet  Union. 

The  overall  objective  of  this  research  program  is  to  systematically  evaluate 
a  wide  variety  of  regional  discriminant  measures,  to  identify  the  reliability  of 
those  discriminants  for  application  to  events  in  the  Soviet  Union,  and  ultimately 
to  formulate  a  winnowing  procedure  utilizing  those  regional  discriminants  which 
can  be  used  to  identify  Soviet  events.  In  prior  investigations  we  have  applied  a 
range  of  regional  signal  measurements,  which  included  new  as  well  as  previ¬ 
ously  proposed  discriminant  measures,  to  the  available  regional  phase  signals 
from  seismic  events  in  tectonic  environments  including  the  Soviet  Union,  the 
western  U.  S.  and  eastern  North  America.  These  measurements  included  time 
domain  phase  amplitude  comparisons  and  spectral  differences  for  events  in 
each  of  the  tectonic  environments.  The  resuits  of  these  studies  have  sug¬ 
gested  that  a  number  of  such  regional  measurements  may  be  useful  in  distin¬ 
guishing  undergound  nuclear  explosions  from  other  seismic  events  including 
earthquakes  and  non-nuclear  blasts.  However,  it  is  only  recently  that  suitable 
regional  data  from  a  large  sample  of  Soviet  explosions  and  nearby  regional 
earthquakes  have  become  available  to  permit  testing  of  those  procedures  on  a 
relevant  database. 

Therefore,  during  this  first  year,  we  have  directed  efforts  at  expanding  the 
database  of  events  from  the  Soviet  Union  and  comparable  events  from  nearby 
border  regions.  This  database  currently  includes  28  East  Kazakh  underground 
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nuclear  explosions  and  32  regional  earthquakes  recorded  at  CD3N  stations 
(primarily  at  station  WMQ).  In  addition,  the  database  of  Soviet  events  includes 
nine  East  Kazakh  explosions  recorded  at  the  Soviet  IRIS  stations,  and  we  are 
in  the  process  of  recovering  a  comparable  sample  of  regional  earthquakes 
recorded  at  these  same  stations.  Additional  regional  recordings  of  Soviet 
underground  nuclear  explosions  obtained  from  these  stations  are  also  being 
incorporated  into  the  database  as  they  become  available. 

We  have  performed  several  analyses  on  these  regional  recordings  con¬ 
centrating  primarily  on  determining  the  detection  capability  of  the  stations  and 
the  amplitude  and  spectral  characteristics  of  the  Lg  and  regional  P  signals. 
Comparisons  of  Lg  versus  regional  P  signal  levels  have  been  analyzed  for 
differences  related  to  source  type.  The  latter  included  investigation  of  fre¬ 
quency  dependence  in  the  relative  excitation  of  Lg  and  regional  P  from  the 
explosion  and  earthquake  sources.  We  also  studied  differences  in  the  relative 
excitation  of  Lg  and  regional  P  signals  between  the  two  principal  explosion  test 
areas  at  East  Kazakh:  Shagan  River  and  Degelen  Mountain. 

Finally,  we  made  some  preliminary  investigations  of  the  effects  of  propa¬ 
gation  path  differences  on  regional  signals  which  could  eventually  be  useful  for 
making  appropriate  adjustments  to  regional  phase  measurements  in  discrimina¬ 
tion  analyses.  In  particular,  since  co-located  earthquake  and  explosion 
sources  do  not  always  exist,  propagation  paths  to  regional  monitoring  stations 
from  explosions  and  comparable  earthquake  sources  may  be  different.  To  get 
at  source-dependent  excitation  differences  more  directly,  signal  measurements 
should  be  compensated  for  attenuation  differences.  In  the  current  study  we 
performed  analyses  to  determine  the  extent  to  which  Lg  signal  differences  at 
the  available  CDSN  and  Soviet  IRIS  stations  could  be  explained  by  attenuation 
differences  for  the  propagation  paths  to  the  stations  from  common  events.  The 
Lg  signal  differences  from  such  common  events  were  used  to  derive  effective 
Q  values  representative  of  the  various  propagation  paths. 

6.2  Conclusions 

The  studies  descnbed  in  this  report  have  revealed  many  interesting 
charactenstics  of  the  behavior  of  regional  signals  from  Soviet  explosions  and 
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comparable  earthquakes,  in  particular,  the  larger  database  of  events  recorded 
at  station  WMQ,  which  was  analyzed  in  this  study,  again  showed  that  Lg/P 
ratios  for  regional  earthquakes  were  relatively  enriched  at  high  frequencies 
compared  to  similar  explosions.  These  Lg/P  ratios  were  found  to  be  intermin¬ 
gled  at  frequencies  near  1  Hz  but  completely  separated  above  2  Hz  for  SR 
explosions  and  regional  earthquakes.  This  observation  continues  to  have  the 
potential  for  becoming  a  valuable  discriminant,  but  effects  of  propagation  path 
differences  between  events  with  different  source  types  need  further  study. 

Another  interesting  conclusion,  based  on  analyses  of  WMQ  recordings  of 
SR  and  DM  explosions,  was  that  the  regional  P  signals  from  SR  explosions  are 
relatively  stronger  than  from  DM  explosions  with  similar  Lg  signals.  This 
appears  to  be  opposite  of  Nuttli's  finding  based  on  comparison  of  teleseismic  P 
from  the  two  source  areas.  Spectral  analyses  of  the  WMQ  signals  from  SR 
and  DM  explosions  indicated  a  tendency  for  the  DM  explosions  to  have  some¬ 
what  larger  Lg/P  ratios  at  higher  frequencies,  which  would  make  them  appear 
more  earthquake-like.  However,  this  effect  was  relatively  small  compared  to 
the  differences  between  explosion  and  earthquake  spectral  ratios  determined  in 
the  discrimination  studies.  Therefore,  it  does  not  severely  limit  the  discrimina¬ 
tion  potential  of  such  measurements. 

An  issue  which  will  require  considerable  study  before  discrimination  results 
can  be  extrapolated  into  uncalibrated  areas  is  the  effect  of  propagation  path 
differences  and  attenuation  on  the  regional  discriminant  measurements.  It  has 
been  argued  that  this  might  be  a  contributing  factor  to  the  Lg/P  ratio  differences 
between  explosions  and  earthquakes  observed  at  WMQ,  which  were  described 
above.  To  better  understand  these  effects,  we  derived  Q  values  representative 
of  Lg  attenuation  for  the  path  from  East  Kazakh  to  WMQ  and  the  other  regional 
stations.  For  the  path  to  WMQ  we  found  a  Q0  value  of  452  representing 
effective  Lg  attenuation.  Q0  values  to  the  other  stations  derived  from  Lg  magni¬ 
tude  residuals  for  common  events  ranged  from  a  low  value  of  428  for  IRIS  sta¬ 
tion  GAR  to  a  high  value  of  761  for  IRIS  station  OBN.  It  is  envisioned  that 
observations  like  these  may  eventually  be  useful  in  adjusting  discriminant 
measurements  for  propagation  differences. 
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In  addition,  some  other  observations  made  during  the  course  of  this  study 
seem  noteworthy.  We  found,  in  general,  that  East  Kazakh  explosions  and 
regional  earthquakes  frequently  produced  strong  Lg  signals  which  were 
recorded  at  CDSN  and  Soviet  IRIS  stations  out  to  ranges  in  excess  of  2000 
km.  For  the  triggered  CDSN  stations,  the  most  complete  sample  of  data  was 
available  from  the  station  WMQ  nearest  the  East  Kazakh  test  site  which  has 
recorded  explosion  signals  to  magnitudes  as  low  as  4.8  mb.  P-wave  signals  at 
WMQ  from  East  Kazakh  explosions  were  seen  to  be  quite  complex.  Spectral 
analyses  of  these  regional  P  signals  revealed  that  selection  of  the  window 
length  appeared  to  produce  some  small  effects  on  computed  spectra,  primarily 
at  lower  frequencies.  Long  windows  encompassing  all  P  energy  produced 
spectra  which  were  somewhat  larger  than  the  corresponding  spectral  levels 
obtained  for  short  windows  including  only  the  initial  P.  As  a  result,  Lg/P  spec¬ 
tral  ratios  obtained  using  the  short  windows  tended  to  be  lower  and  somewhat 
more  earthquake-like  than  ratios  based  on  the  longer  P-wave  windows.  Possi¬ 
ble  causes  for  this  observation  need  additional  study.  CDSN  station  HIA  at  a 
range  of  more  than  2900  km  from  the  test  site  also  recorded  strong  Lg  signals 
from  many  of  the  larger  East  Kazakh  explosions.  However,  the  useful  fre¬ 
quency  band  of  these  signals  was  found  to  be  limited  to  rather  low  frequencies. 
Unfortunately,  other  CDSN  stations  have  not  recorded  useful  Lg  signals  from 
Soviet  explosions  apparently  because  their  normal  trigger  level  is  set  higher 
than  the  amplitudes  generated  by  the  signals  at  these  ranges. 

Much  of  our  analysis  of  the  Soviet  IRIS  station  data  is  preliminary  at  this 
time  since  we  are  still  developing  the  earthquake  database.  However,  we  have 
found  that  stations  ARU  and  GAR  frequently  record  strong  Lg  and  regional  P 
signals  from  East  Kazakh  explosions  down  to  fairly  low  magnitudes.  Spectral 
analysis  of  these  signals  indicates  that  the  maximum  signal-to-noise  level  for 
these  events  is  in  a  oassband  near  I  Hz.  In  fact,  we  found  that,  by  band-pass 
filtering  the  ARU  records  in  a  frequency  band  from  0.8  to  1.6  Hz,  Lg  signal  lev¬ 
els  could  be  enhanced  to  three  times  the  noise  level  for  a  SR  explosion  with  a 
magnitude  of  only  about  3.8  mD.  Similar  enhancement  also  is  attainable  for  Lg 
signals  recorded  at  IRIS  station  GAR.  However,  at  the  more  distant  Soviet 
IRIS  stations,  the  Lg  signals  from  East  Kazakh  explosions  have  spectral  peaks 
at  frequencies  less  than  1  Hz  and  drop  off  rapidly  to  the  noise  level  at  higher 
frequencies.  As  a  result,  it  appears  ttu  t  stations  considerably  closer  than  O BN 
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and  KIV  will  probably  be  required  to  obtain  Lg  signals  at  frequencies  above  1 
Hz.  This  may  impede  some  types  of  regional  discriminant  measures  which 
require  broadband  spectral  estimates  of  signal  strength. 
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Seixmographic  Station 
University  of  California 
Berkeley,  CA  94720 


Prof.  Alan  Kafka 

Department  of  Geology  &  Geophysics 

Boston  College 

Chestnut  Hill,  M A  02167 


Dr.  Richard  LaCoss 
MIT-Eincoln  Laboratory 
M-200B 
P.  O.  Box  73 

Lexington,  MA  02173-0073  (3  copies) 
Prof  bred  K.  Lamb 

University  of  Illinois  at  Urbana-Champaign 
Department  of  Physics 
1110  West  Green  Street 
Urbana,  IL  61 H01 

Prof.  Charles  A.  Langston 
Geosciences  Department 
403  Deike  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16X02 


Prof.  Thorne  La}' 

Institute  of  Tectonics 

Earth  Science  Board 

University  of  California,  Santa  Cruz 

Santa  Cruz,  CA  95064 

Prof.  Arthur  Lemer-Lam 
Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


Dr.  Christopher  Lynnes 
Teledyne  Geotech 
314  Montgomery  Street 
Alexandria,  VA  22314 


Prof.  Peter  Mai  in 

University  of  California  at  Santa  Barbara 
Institute  for  Crustal  Studies 
Santa  Barbara,  CA  9310ft 


Dr.  Randolph  Martin,  III 
New  England  Research,  Inc. 

76  Olcott  Drive 

White  River  Junction,  VT  05001 


Dr.  Gary  McCanor 
Mission  Research  Corporation 
735  Slate  Street 
P.O.  Drawer  719 

Santa  Barbara,  CA  93102  (2  copies) 

Prof.  Thomas  V.  McEvilly 
Seismographic  Station 
University  of  California 
Berkeley,  CA  94720 


Dr.  Keith  L.  McLaughlin 
S-CUBED 

A  Division  of  Maxwell  Laboratory 
P.O.  Box  1620 
La  Jolla,  CA  92038-1620 

lYof.  William  Menke 
Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


Stephen  Miller 
SRI  International 
333  Ravenswood  Avenue 
Box  AF  1 1 6 

Menlo  Park,  CA  94025-3493 


Prof.  Bernard  Minster 
IGPP.  A -02 5 

Senpps  Institute  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92023 

Prof.  Brian  J.  Mitchell 

Department  of  Earth  &  Atmospheric  Sciences 
St.  1-ouis  University 
St.  Louis,  MO  63156 


Mr.  Jack  Murphy 

S-CUBED,  A  Division  of  Maxwell  Laboratory 
1 1  S(X)  Sunrise.  Valiev  Drive 
Suite  1212 

Keston,  VA  22091  (2  copies) 

Dr.  Bao  Nguyen 
GL/LWII 

Ilanscom  AFB,  MA  01731-5000 


Prof.  John  A.  Orcutt 
IGPP,  A-025 

Scripps  Institute  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92093 

Prof.  Keith  Priestley 
University  of  Cambridge 
Bullard  Labs,  Dept,  of  Earth  Sciences 
Madinglcy  Rise,  Madingley  Rd. 
Cambridge  CB3  OEZ,  ENGLAND 

Prof.  Paul  G.  Richards 
Lamont  Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


Dr.  Wilmer  Rivers 
Teledyne  Geotcch 
314  Montgomery  Street 
Alexandria,  VA  22314 


Prof.  Charles  G.  Sammis 
Center  for  Earth  Sciences 
University  of  Southern  California 
University  Park 
Los  Angeles,  CA  90089-0741 

Prof.  Christopher  H.  Scholz 
Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Pali^dcs,  NY  10964 
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Thomas  J.  Sereno,  Jr. 

Science  Application  Inl’l  Corp. 
10260  Campus  Point  Drive 
San  Diego,  CA  92121 


Prof.  David  G.  Simpson 
Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


Dr.  Jeffrey  Stevens 
S-CUBED 

A  Division  of  Maxwell  Laboratory 
P.O.  Box  1620 
La  Jolla,  CA  92038-1620 

Prof.  Brian  Stump 

Institute  for  the  Study  of  Earth  &.  Man 
Geophysical  Laboratory 
Southern  Methodist  University 
Dallas,  TX  75275 

Prof.  Jeremiah  Sullivan 

University  of  Illinois  at  Urbana-Champaign 

Department  of  Physics 

1 1 10  West  Green  Street 

Urbana,  IL  61801 

Prof.  Clifford  Thurber 
University  of  Wisconsin-Madison 
Department  of  Geology  &  Geophysics 
1215  West  Dayton  Street 
Madison,  WS  53706 

Prof.  M.  Nafi  Toksoz 
Earth  Resources  Lab 
Massachusetts  Institute  of  Technology 
42  Carlcton  Street 
Cambridge,  MA  02142 

Prof.  John  E.  Vidale 

University  of  California  at  Santa  Cruz 

Seismological  Laboratory 

Santa  Cruz,  CA  95064 


Prof.  Terry  C.  Wallace 
Department  of  Geosciences 
Building  #77 
University  of  Arizona 
Tucson,  AZ  85721 

Dr.  Raymond  Willeman 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 

Dr.  Lorraine  W^i f 
GL/LWH 

Hanscom  Al’B,  MA  01731-5000 


Dr.  Monem  Abdel-Gawad 
Rockwell  International  Science  Center 
1049  Camino  Dos  Rios 
Thousand  Oaks,  CA  91360 


l’rof.  Kcnti  Aki 

Center  for  Earth  Sciences 

University  of  Southern  California 

University  Park 

Los  Angeles,  CA  90089-0741 

Prof.  Shelton  S.  Alexander 
Geosciences  Department 
403  IXike  Building 
'Hie  Pennsylvania  State  University 
University  Park,  PA  16802 

Dr.  Kenneth  Anderson 
BBNSTC 
Mail  Stop  14/1 B 
Cambridge,  MA  02238 


Dr.  G.A.  Bollinger 
Department  of  Geological  Sciences 
Virginia  Polytechnical  Institute 
21044  Dcrnng  1  fall 
Blacksburg,  VA  24061 

Dr.  Stephen  Bratt 
Center  for  Seismic  Studies 
1300  North  1 7th  Street 
Suite  1450 

Arlington,  VA  22209 

Michael  Browne 
Teledyne  Geotech 
3401  Shiloh  Road 
Garland,  TX  75041 


Mr.  Roy  Burger 
1221  Serry  Road 
Schenectady,  NY  12309 


Dr.  Ralph  Archuleta 
Department  of  Geological  Sciences 
University  of  California  at  Santa  Barbara 
Santa  Barbara,  CA  93102 


Dr.  Jeff  Barker 

Department  of  Geological  Sciences 
State  University  of  New  York 
at  Binghamton 
Vestal,  NY  13901 

Dr.  Susan  Beck 

Department  of  Geosciences,  Bldg  #  77 
University  of  Arizona 
Tucson,  AZ  85721 


Dr.  T.J.  Bennett 
S-CUBED 

A  Division  of  Maxwell  Laboratory 
1 1800  Sunrise  Valley  Drive,  Suite  1212 
Reston,  VA  22091 

Mr.  William  J.  Best 
907  Westwood  Drive 
Vienna,  VA  22180 


Dr.  N.  Biswas 
Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 


Dr.  Robert  Bumdge 
Schlumberger-Doll  Research  Center 
Old  Quarry  Road 
Ridgefield,  CT  tk>877 


Dr.  W.  Winston  Chan 
Tcledyne  Geotech 
314  Montgomery  Street 
Alexandria,  VA  22314-1581 


Dr.  Theodore  Cherry 
Science  Horizons,  Inc. 

710  Encinitas  Blvd.,  Suite  200 
Encinitas,  CA  92024  (2  copies) 


Prof.  Jon  F.  Claerbout 
Department  of  Geophysics 
Stanford  University 
Stanford,  CA  94305 


Prof.  Robert  W.  Clayton 

Seismologicai  Laboratory 

Division  of  Geological  Sc  Planetary  Sciences 

California  Institute  of  Technology 

Pasadena,  CA  91125 

Prof.  F.  A.  Dahlen 

Geological  and  Geophysical  Sciences 

Princeton  University 

Princeton,  NJ  08544-0636 
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Prof.  Adam  Dziewonski 
Hoffman  Laboratory 
1  larvard  University 
20  Oxford  St 
Cambridge,  MA  02138 

.  Prof.  John  Ebel 

Department  of  Geology  &  Geophysics 
Boston  College 
Chestnut  Mill,  MA  02167 

« 


Eric  Fielding 
SNEE  Midi 

IMSTOP 

Corned  Umvci'Mity 
Ithaca,  NY  14853 

Prof.  Donald  Forsyth 
Department  of  Geological  Sciences 
Brown  University 
Providence,  RI  02912 


Dr.  Cliff  Frolich 
Institute  of  Geophysics 
8701  North  Mopac 
Austin,  IX  78759 


Dr.  Anthony  Gangi 
Texas  A&M  University 
Department  of  Geophysics 
College  Station,  TX  77843 


Dr.  Freeman  Gilbert 
IGPP,  A -02 5 

Scripps  Institute  of  Oceanography 
University  of  California 
La  Jolla,  CA  92093 

Mr.  Edward  Giller 
Pacific  Sierra  Research  Corp. 

1401  Wilson  Boulevard 
Arlington,  VA  22209 


Dr.  Jeffrey  W.  Given 
SA1C 

.  10260  Campus  Point  Drive 
San  Diego,  CA  92121 

Prof.  Roy  Greenfield 
Geosciences  Department 
403  Dcike  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Dan  N.  I  lagedorn 
Batlelle 

Pacific  Northwest  Laboratories 
Battellc  Boulevard 
Richland,  WA  99352 

Kevin  Hutchenson 
Department  of  Larth  Sciences 
St.  Louis  University 
3507  Laclede 
St.  Louis,  MO  63103 

Dr.  Hans  Israelsson 
Center  for  Seismic  Studies 
1300  N.  1 7th  Street,  Suite  1450 
Arlington,  VA  22209-2308 


Prof.  Thomas  II.  Jordan 
Department  of  Earth,  Atmospheric 
and  Planetary  Sciences 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 

Robert  C.  Kemerait 
ENSCO,  Inc. 

445  Pineda  Court 
Melbourne,  EL  32940 


William  Kikendall 
Teledyne  Geotech 
3401  Shiloh  Road 
Garland,  TX  75041 


Prof.  Leon  Knopoff 

University  of  California 

Institute  of  Geophysics  &  Planetary  Physics 

Los  Angeles,  CA  90024 

Prof.  L.  Timothy  Long 
School  of  Geophysical  Sciences 
Georgia  Institute  of  Technology 
Atlanta,  GA  30332 


Prof.  Art  McGarr 
Mail  Stop  977 
Geological  Survey 
345  Middieficld  Rd. 
Menlo  Park,  CA  94025 

Dr.  George  Mellman 
Sierra  Geophysics 
11255  Kiikland  Way 
Kirkland,  WA  98033 


Prof  John  Nabelck 
College  of  Oceanography 
Oregon  State  University 
Corvallis,  OR  97331 


Dr.  Susan  Schwartz 
Institute  of  Tectonics 
1156  High  St. 

Santa  Cruz,  CA  95064 


Prof.  Geza  Nagy 

University  of  California,  San  Diego 
Department  of  Ames,  MS.  13-010 
La  Jolla,  CA  92093 


Prof.  Amos  Nur 
Department  of  Geophysics 
Stanford  University 
Stanford,  CA  94305 


Prof.  Jack  Oliver 
Department  of  Geology 
Cornell  University 
Ithaca,  NY  14X50 


Prof.  Robert  Phinncy 
Geological  &  Geophysical  Sciences 
Princeton  University 
Princeton,  NJ  0X544-0636 


Dr.  Paul  Pomeroy 
Rondout  Associates 
P.O.  Box  224 
Stone  Ridge,  NY  124X4 


Dr.  Jay  Pulli 
RADIX  System,  Inc. 

2  Taft  Court,  Suite  203 
Rockville,  MD  20X50 


Dr.  Norton  Rimer 
S-CUBED 

A  Division  of  Maxwell  Laboratory 
P.O.  Box  1620 
La  Jolla,  C A  9203X-1620 

Prof.  Larry  J.  Ruff 
Department  of  Geological  Sciences 
1006  C.C.  Little  Building 
University  of  Michigan 
Ann  Arbor,  MI  4X109-1063 

Dr.  Richard  Sailor 
TASC  Inc. 

55  Walkers  Brook  Drive 
Reading,  MA  01X67 


John  Sherwin 
Teledvne  Geotech 
3401  Shiloh  Road 
Garland,  TX  75041 


Dr.  Matthew  Sibol 
Virginia  Tech 

Seismological  Observatory 
4044  Derring  Hall 
Blacksburg,  VA  24061-0420 

Prof.  Robert  Smith 
Department  of  Geophysics 
University  of  Utah 
1400  East  2nd  South 
Salt  Lake  City,  UT  X41 12 

Dr.  Stewart  W.  Smith 
Geophysics  AK-50 
University  of  Washington 
Seattle,  WA  98195 


Dr.  George  Sutton 
Rondout  Associates 
P.O.  Box  224 
Stone  Ridge,  NY  12484 


Prof.  L.  Sykes 

Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


Prof.  Pradeep  Talwani 
Department  of  Geological  Sciences 
University  of  South  Carolina 
Columbia,  SC  29208 


Prof.  Ta-liang  Teng 
Center  for  Earth  Sciences 
University  of  Southern  California 
University  Park 
Los  Angeles,  CA  90089-0741 

Dr.  R.B.  Tittmann 

Rockwell  International  Science  Center 

1049  Camino  Dos  Rios 

P.O.  Box  1085 

Thousand  Oaks,  CA  91360 


Dr.  Gregory  van  der  Vink 
IRIS,  Inc. 

1016  North  Fort  Myer  Drive 
Suite  1440 

Arlington,  VA  22209 

Professor  Daniel  Walker 
University  of  lawaii 
Institute  of  Geophysics 
Honolulu,  HI  96822 


William  R.  Walter 
Seismological  Laboratory 
University  of  Nevada 
Reno,  NV  89557 


Dr.  Gregory  Wojcik 
Weidlinger  Associates 
4410  El  Camino  Real 
Suite  1 10 

Los  Altos,  CA  94022 

Prof.  John  H.  Woodhousc 
Hoffman  Laboratory 
Harvard  University 
20  Oxford  St. 

Cambridge,  MA  02138 

Prof.  Francis  T.  Wu 
Department  of  Geological  Sciences 
State  University  of  New  York 
at  Binghamton 
Vestal,  NY  13901 

Dr.  Gregory  B.  Young 
ENSCO,  Inc. 

5400  Port  Royal  Road 
Springfield,  VA  22151-2388 
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Dr  Ralph  Alewine  111 
DARPA/NMRO 
1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 


Mr.  James  C.  Battis 
Gl/LWH 

1  lanscom  AFB,  MA  01731-5000 


Dr.  Robert  Blandford 
DAKPA/NMRO 
1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 


Brie  Chael 
Division  9241 
Sandia  Laboratory 
Albuquerque,  NM  K71H5 


Dr.  T.  Hanks 
USGS 

Nat'l  Earthquake  Research  Center 
345  Middleficld  Road 
Menlo  Park,  CA  94025 

Dr.  James  1  lannon 

Lawrence  Livermore  Nat'l  Laboratory 
P.O.  Box  808 
Livermore,  CA  94550 


Paul  Johnson 

ESS-4,  Mail  Stop  J979 

Los  Alamos  National  Laboratory 

Los  Alamos,  NM  87545 


Janet  Johnston 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Dr.  John  J.  Cipar 

GL/LWH 

Hanscom  ABB,  MA  01731-5000 

Dr.  Katharine  Kadinsky-Cade 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 

Cecil  Davis 

Group  P  15,  Mail  Stop  D406 

P.O.  Box  1663 

Los  Alamos  National  Laboratory 

Los  Alamos,  NM  87544 

Ms.  Ann  Kerr 

IGPP,  A-025 

Scripps  Institute  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92093 

Mr.  JelT  Duncan 

Office  of  Congressman  Markey 

2133  Rayburn  House  Bldg. 

Washington,  DC  20515 

Dr.  Max  Koontz 

US  Dept  of  Encrgy/DP  5 

Forrcslal  Building 

1000  Independence  Avenue 
Washington,  DC  20585 

Dr.  Jack  Bvemden 

USGS  -  Earthquake  Studies 

345  Middleficld  Road 

Menlo  Park,  CA  94025 

Dr.  W.H.K.  Lee 

Office  of  Earthquakes,  Volcanoes, 

&  Engineering 

345  Middle  fie  Id  Road 

Menlo  Park,  CA  94025 

Art  I'Yankel 

USGS 

922  National  Center 

Reston,  VA  22092 

Dr.  William  Leith 

U.S.  Geological  Survey 

Mail  Stop  928 

Reston,  VA  22092 

Dr.  Richard  Lewis 

Director,  Earthquake  Engineering  &  Geophysics 
U.S.  Army  Corps  of  Engineers 
Box  63 1 

Vicksburg,  MS  39180 


Dr.  Dale  Glover 
DIA/DT- 1 B 

Washington,  DC  20301 


fl 


James  l'.  Lewkowicz 
GL/LWI  1 

1  Ianscom  AFB,  MA  01731-5000 


»  Mr.  Alfred  Liebennan 

ACDA/VI-OA'State  Department  Bldg 
Room  5726 

r  320  -  21  st  Street,  NW 
Washington,  DC  20451 

Stephen  Mangino 
GL/LWH 

1  Ianscom  AFB,  MA  01731-5000 


Mr.  Chris  Paine 
Office  ol  Senator  Kennedy 
SR  315 

United  States  Senate 
Washington,  DC  20510 

Colonel  Jerry  J.  Perrizo 
AFOSR/NP,  Building  410 
Bolling  AFB 

Washington,  DC  20332-6448 


Dr.  Frank  F.  Pilotte 

HQ  AFTAGTT 

Patrick  AFB,  FL  32925-600! 


Dr.  Robert  Masse 
Box  25046,  Mail  Stop  967 
Denver  Federal  Center 
Denver,  CO  80225 


Katie  Poley 
C1A-OSWR/NED 
Washington,  DC  20505 


Art  McGarr 

U.S.  Geological  Survey,  MS-977 
345  Middle-field  Road 
Menlo  Park,  CA  94025 


Richard  Morrow 
ACDA/VI,  Room  5741 
320  21st  Street  N.W 
Washington,  DC  20451 


Dr.  Keith  K.  Nakanishi 
Lawrence  Livermore  National  Laboratory 
P.O.  Box  808,  L-205 
Livermore,  CA  94550 


Dr.  Carl  Newton 

Los  Alamos  National  Laboratory 

P.O.  Box  1663 

Mail  Stop  C335,  Group  ESS-3 
Los  Alamos,  NM  87545 
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Dr.  Kenneth  H.  Olsen 
Los  Alamos  Scientific  Laboratory 
-  P.O.  Box  1663 
Mail  Stop  D-406 
Los  Alamos,  NM  87545 

Howard  J.  Patton 

Lawrence  Livermore  National  Laboratory 
P.O.  Box  808,  L-205 
Livermore,  CA  94550 


Mr.  Jack  Rachlin 
U.S.  Geological  Survey 
Geology,  Rm  3  Cl 36 
Mail  Stop  928  National  Center 
Reston,  VA  22092 

Dr.  Robert  Reinke 
WL/NTESG 

Kirtland  AFB,  NM  87117-6008 


Dr.  Byron  Ristvet 

HQ  DNA,  Nevada  Operations  Office 

Attn:  NVCG 

P.O.  Box  98539 

Las  Vegas,  NV  89193 

Dr.  George  Rothe 
HQ  AFTAC/TTR 
Patrick  AFB,  FL  32925-6001 


Dr.  Alan  S.  Ryall,  Jr. 
DARPA/NMRO 
1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 


Dr.  Michael  Shore 
Defense  Nuclear  Agency/SPSS 
6801  Telegraph  Road 
Alexandria,  VA  22310 


o 


Dr.  Albert  Smith 

Lain  Alamos  National  Laboratory 

L-205 

P.  O.  Box  SOX 
Livermore.  CA  94550 

Donald  L.  Springer 

Lawrence  Livennorc  National  Laboratory 
L-205 

P.  O.  Box  SOS 
Livemiore  ,  CA  94550 

Mr.  Charles  L.  Taylor 
C17LWC 

Hanscom  AFB,  MA  01731-5000 


Dr.  Steven  R.  Taylor 

Lawrence  Livemiore  National  Laboratory 

L-205 

P.  O.  Box  SOS 
Livermore,  CA  94550 

Dr.  Larry  Turnbull 
CIA-OSWR/NLD 
Washington,  DC  20505 


Dr.  Eileen  Vergino 

Lawrence  Livemiore  National  Laboratory 
L-205 

P.  O.  Box  SOS 
Livemiore,  CA  94550 

Dr.  Thomas  Weaver 
Ixs  Alamos  National  Laboratory 
P.O.  Box  i663.  Mail  Stop  C335 
Los  Alamos,  NM  X7545 


J.J  Zucca 

Lawrence  Livermore  National]  Laboratory 
P.  O.  Box  S0K 
Livemiore,  CA  94550 


GUSULL 
Research  Library 

I  lanscom  AFB  ,  MA  01731-5000  (2  copies) 


Secretarv  of  the  Air  Force 
(SAIRD) 

Washington,  DC  20330 


Office  of  the  Secretary  Defense 
DDR  &  E 

Washington,  DC  20330 


IIQDNA 

Attn:  Technical  Library 
Washington,  DC  20305 


D  A  RPA/RMO/R  ETR I E  V  A  L 
1400  Wilson  Boulevard 
Arlington,  VA  22209 


DARPA/RMO/Security  Office 
1400  Wilson  Boulevard 
Arlington,  VA  22209 


Geophysics  Laboratory 
Attn:  XO 

Hanscom  AFB,  MA  01731-5000 


Geophysics  Laboratory 
Attn:  LW 

Hanscom  AFB.  MA  01731-5000 


DARPA/PM 

1400  Wilson  Boulevard 

Arlington,  VA  22209 


Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22314  (5  copies) 


Defense  intelligence  Agency 
Directorate  for  Scientific 
&  Technical  IntelligcnceAttn:  DT1B 
Washington,  DC  20340-6158 


AFTAC/CA 

(ST1NFO) 

Murick  AFR,  FL  32925-600! 

TACTEC  -  nattelle  Memorial  Institute 
00r>  King  Avp, 

Columbus,  011  43201  (KINAL  REPORT  ONLY) 
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CONTRACTOR? 


Dr.  Ramon  Cabre,  S.J. 
Observaturio  S,ui  CaJixlo 

c.isiiia  5iw 

La  Pa/.,  liolivia 


’Prof.  1  lans-Pcter  Jlarjes 
Institute  lor  Geophysik. 
Ruhr  Umversity/Uochum 
•P.O.  Box  10214K 
4(00  Bochum  1 ,  FRG 

Prof,  Eystein  1  luscbyc 
NTNIVNORSAR 
P.O.  Box  .SI 

N-20O7  Kjclicr,  NORWAY 


Prof.  Brian  L.N.  Kennctt 
Research  School  of  Barth  Sciences 
Institute  of  Advanced  Studies 
G.P.O.  Box  4 

Canberra  2001,  AUSTRALIA 

Dr.  Bernard  Massinon 

Societe  Radiomana 

27  rue  Claude  Bernard 

75005  Paris,  FRANCIS  (2  Copies) 


Dr.  Pierre  Mechcler 
Societe  Radiomana 
27  rue  Claude  Bernard 
75005  Paris,  FRANCE 


Dr.  Svein  MykJceltveit 
NTNIVNORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


f 


« 
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FOrtt:  ICN  (Ot  nor:,) 


Dr.  I  ’c ter  Basham 

Earth  Physics  Branch 

Geological  Survey  o!  Canada 

!  Observatory  Crescent 

Ottawa.  C'.aario.  CANADA  K.1A  0Y3 


Dr.  Eekadu  Kebede 
Seismological  Section 
llox  12019 

S-750  Uppsala,  SWEDEN 


Dr  Eduard  Berg 
Institute  o!  Geophysics 
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